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Abstract 



In part I we study quantum modified photon trajectories in a Schwarzschild 
blackhole spacetime. The photon vacuum polarization effect in curved space- 
time leads to birefringence, i.e. the photon velocity becomes c±(5c depending 
on its polarization. This velocity shift then results in modified photon tra- 
jectories. 

We find that photon trajectories are shifted by equal and opposite amounts 
for the two photon polarizations, as expected by the sum rule [H]. Therefore, 
the critical circular orbit at u = 1/3M in Schwarzschild spacetime, is split 
depending on polarization as u = 1/3M±A6{M) (to first order in A), where 
A is a constant found to be ~ 10"'^^ for a solar mass blackhole. Then using 
general quantum modified trajectory equations we find that photons pro- 
jected into the blackhole for a critical impact parameter tend to the critical 
orbit associated with that polarization. We then use an impact parameter 
that is lower than the critical one. In this case the photons tend to the event 
horizon in coordinate time, and according to the affine parameter the pho- 
tons fall into the singularity. This means even with the quantum corrections 
the event horizon behaves in the classic way, as expected from the horizon 
theorem [H]. 

We also construct a quantum modified Schwarzschild metric, which en- 
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compasses the quantum polarization corrections. This is then used to derive 
the photons general quantum modified equations of motion, as before. We 
also show that when this modified metric is used with wave vectors for radi- 
ally projected photons we obtain the classic equations of motion, as expected, 
because radial velocities are not modified by the quantum polarization cor- 
rection. 

In Part II we use the 2+ld Nambu-Jona-Lasino (NJL) model to study the 
superfluid behaviour of two-dimensional quark matter. In previous work, 
[T2] . it was suggested that the high density phase of the 2+ld NJL model 
could be a relativistic gapless thin film BCS superfiuid. In this work we find 
that as we raise the baryon chemical potential (/i) the baryon supercurrent 
jumps from a non-superfluid (zero) phase to a superfluid (non-zero) phase. 
This sharp transition is seen to occur in the region 0.65 < jic < 0.68, which 
was shown in ^2] to be the region of chiral symmetry restoration. In this 
analysis we prove that at high density the 2 + ld NJL model is in a superfluid 
phase. 

We then go on to study the dynamics of the superfluid phase, represented 
by the helicity modulus (T ), which is the constant of proportionality between 
the supercurrent and the gradient of the diquark state function. We find that 
below the temperature associated with lattice size Lt = 4, the system is in a 
non-superfiuid phase, and above Lt = 24 the system is in a superfiuid phase. 
We also find a possible 2nd order transition at Lt ^ 6, which corresponds 
to the critical point as described by Kosterlitz and Thouless' theory of 2D 
critical systems with U{1) global symmetry - such as the XY model. 
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Preface 



In this work we will be dealing with quantum particles under extreme condi- 
tions. This will include the study of quantum modified photon trajectories in 
a highly curved background spacetime (Schwarzschild blackhole spacetime) 
and high density quark matter in two dimensional films. 

Part 1: Quantum Modified Null Trajectories in Schwarzschild Space- 
time 

In this part we will be studying the quantum modifications to null geodesies 
in Schwarzschild spacetime, in particular the modifications to the critical 
stable orbit. In Chapter^ we will give an introduction to the the vacuum 
polarization efi"ect that leads to the modification of photon trajectories in 
general relativity. In Chapter [21 we will give an introduction to the derivation 
of the orbit equation in Schwarzschild spacetime, including the techniques 
involved in solving it and the interpretation of results, i.e. for critical stable 
orbits and trajectories to the singularity. In Chapter El we will then provide 
an introduction to the techniques of geometric optics, i.e. the derivation of 
the light-cone and geodesic equation from Maxwell's equations of motion in 
curved spacetime (with and without the quantum corrections). We will also 
give a brief overview of previous theories associated with quantum modified 
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null geodesies, such as the horizon theorem and the polarization sum rule. In 
Chapter m we will go onto using theses techniques in order to determine the 
quantum modified critical circular orbits, and the changes in the associated 
critical impact parameters. These finding will then be extended, by deriving a 
general quantum modified equation of motion. Using this, with the modified 
impact parameters, we will study the modifications to general null geodesies 
as they tend to the critical stable orbits. Also, using the quantum modified 
general equation of motion, we will show that a general geodesic aimed into 
the black hole behaves in the classical way around the event horizon, as 
described by the horizon theorem. In Chapter El we will encompass these 
results into a quantum modified Schwarzschild metric. 

Part 2: Study of the 2+ld NJL Gapless Superfiuid 

In this part we will be studying high density quark matter in two dimensional 
films using the 2+ld Nambu-Jona-Lasino (NJL) model, in particular to iso- 
late a critical gapless BCS superfiuid phase, which was suspected to exist in 
[T2] . In Chapter Owe will set the scene for this work by giving an overview 
of quantum chromodynamics, its symmetries and aspects of its phase dia- 
gram. In Chapter IHl we will give a brief introduction to the NJL model; then 
we will go on to discuss the techniques and results of the previous simula- 
tions of the NJL model, given in 2+ld [12], i.e. simulation of the diquark 
condensate through the introduction of diquark sources, and the evidence of 
chiral symmetry restoration and the non-zero baryon density. In Chapter [HI 
we will extend the baryon density to a baryon three current by the use of 
Ward identities, these are then implemented into the simulation. We also 
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introduce a spatially varying diquark source, referred to as a twisted source. 
In this way a gradient in the diquark pair wave function is introduced, which 
forces a flow of the baryon current, which will then be measured. We will 
then explore the behaviour of the 3-current (using the helicity modules, T, 
Sec. 19. 2j) with variations in spatial volume, temperature, and variations in 
the diquark source, in order to isolate the superfluid (and non-superfluid) 
state of the 2+ld NJL model. In Chapter [HI we will study the variation 
of the helicity modulus with temperature; which will be done in order to 
determine the critical point of the system, i.e. the point where vortex and 
antivortex pairs come together to form the superfluid phase (as predicted in 
the condensed matter study of 2 dimensional systems, e.g. the XY model). 
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Trajectories in Schwarzschild 
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Chapter 1 
Introduction 

1.1 General Relativity 

Since the birth of special relativity, in 1905, the nature of space and time has 
been demoted to a relative entity, known as spacetime, which is stretched and 
contracted depending on an observer's frame of reference, while the speed of 
light, c, has taken the pedestal of an absolute and universal speed limit, 
unaffected by any transformation of reference frame. From this emerged a 
generalized theory of relativity, which portrayed the gravitational field in a 
new and revolutionary way: where it didn't depend on a propagating field 
but on the nature of spacetime itself. In this view matter (or energy) is 
said to curve and modify the surrounding spacetime, this then results in 
photons and particles tracing out shortest paths between two points, known 
as geodesies. Therefore, gravitational forces become a manifestation of the 
curved spacetime due to the presence of matter |U[S]. In this general rela- 
tivistic framework spacetime is described by the metric g^i, and the motion 
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CHAPTER 1. INTRODUCTION 3 

of particles are described by the interval equation: 

> Time-like (c < 1) 
k"^ = g^^k^^k" ^ =0 Light-like (c = 0) (1.1) 

< Space-like (c > 1) 

where = For flat spacetime or a local inertia! frame (LIF), where 

g^u is replaced by the diagonal Minkowski metric rj^^ = (1,-1,— 1,-1), the 
interval equation becomes: 

,9 dt'^ dx'^ dy"^ dz"^ ,^ 

*= = ri,.k-e = (1.2) 

for c = 1. Apart from resolving the problems associated with Newtonian me- 
chanics, such as describing the perihelion advance of Mercury, the strongest 
aspect of general relativity was its predictive power. One of its most radical 
claims, and the building blocks of the theory itself, was that gravitational 
fields affect radiation, which was then confirmed through the observation 
of starlight deflection by the sun. From this emerged some profound and 
fantastic possibilities such as black holes and gravitational lensing. The 
gravitational eff'ect on light rays also leads to the possibility that photons 
could follow stable orbits around stars (discussed in chapter E]), and it's this 
possibility in which we will be interested. 

1.2 QED in a Curved Spacetime 

Even though photon trajectories are modified in a curved spacetime, and the 
resulting curved paths are described by general relativity, this bending of 
light was, for a long time, considered to have no effect on the velocity of the 

photon. This view shifted slightly when, in 1980, Drummond and Hathrell[T] 

^For photons this becomes = written in terms of the affine parameter A 
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Figure 1.1: First order, a, vacuum polarization Feynman diagram contribut- 
ing to the photon propagator 

proposed that a photon propagating in a curved spacetime may, depending on 
its direction and polarization, travel with a velocity that exceeds the normal 
speed of light c. This change in velocity would then result in trajectories 
other then the ones described by "classical" general relativity. This effect is 
simply described as a modification of the light cone in a LIF: 

= Vf..k''k' = ^ (r/^, + (1.3) 

Where a is the fine structure constant and a^iy{R) is a modification to the 
metric that depends on the Riemann curvature at the origin of the LIF. This 
correction is seen to arise from photon vacuum polarizations in a curved 
space time. Fig. 11.11 Qualitatively it can be thought of as a photon splitting 
into a virtual e~^e~ pair, so at the quantum level it is characterized by the 
Compton wavelength Ac; then, when this quantum cloud of size 0{Xc) passes 
through a curved spacetime its motion would be affected differently to that 
described by general relativity, possibly in a polarization-dependent way [H]. 

This effect of vacuum polarization is considered through the effective 
action: 

S = Si + S2 (1.4) 
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where 5*1 is the Maxwell electromagnetic action in curved space time: 

F^u = d^A,~d,A^ (1.5) 

and 5*2 is the part of the action that incorporates the effects of virtual electron 
loops in a background gravitational field. As we are only concerned with 
the propagation of individual photons it must be quadratic in A^, and the 
constraint of gauge invariance then implies that it must depend on F^y rather 
than Afj_. Also, as the virtual loops give the photon a size of Ac, 5*2 can be 
expanded in powers of = m~^, thus the lowest term in the expansion 
would be of order m"^, which is the term corresponding to one electron loop. 
With these constraints there are only four independent gauge invariant terms, 
which can be chosen to be: 

S2 = ^Jd^xV^iaRF^yF^^'^ + bR^yF^^'^F^ 

+ cR^y^rF'^'F''^ + dD^F''''D^F^) (1.6) 

The first three terms represent a direct coupling of the electromagnetic field 
to the curvature, and they vanish in flat-spacetime. The fourth, however, 
is also applicable in the case of flat-spacetime, and represents off-mass-shell 
effects in the vacuum polarization. In the values for a, 6, c, and d have 
been determined to 0{e^). The constant d is obtained by comparing the 
coefficient of the renormalized flat-spacetime photon propagator associated 

with the Feynman diagram in Fig. II. if to the result of the same order given 

^The photon propogator with the vacuum polarization in flat-spacetime is given (in the 
Feynman gauge) by: ^ ^ ^ + ^I^^., where /^^ = {Vf^-^q^ - qi,qu){l - ef^;^ + • ■ •) 
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by the effective action S; and a, b, and c are obtained by comparing the 
coefficients of the coupling of a graviton to two photons'^ to the same result 
obtained from S2. In this way the constants are given as: 

5 a 26 a 2 a 24 a 

a = = c= d = (1.7) 

720 TT 720 TT 720 vr 720 tt ^ ^ 

Then, as the equations of motion for the electromagnetic field are given by: 

(1.8) 



using the modified action ()1.4p we find: 

D,F'^^ + ^ = (1.9) 

From this we can see that D^j^F^^ is of C(e^), therefore, the term with coef- 
ficient d in Eqn. ()1.6|) will be of O(e^), hence we can omit it from the final 
equation of motion. In this way we find [2|: 

D^F^^ - ^[2bR^xD^F^'' + Acg^^R^rXpD^F^P] = (1.10) 

which is the Maxwell equation in curved spacetime, incorporating the cou- 
pling of curvature with vacuum polarization effects. Using this modified 
Maxwell equation and the methods of geometric optics (described in Chap- 
ter E]) it is possible to derive the quantum modified light cone and geodesic 
equations. Then, using these, we are able to determine the quantum modified 

trajectories in curved spacetime. 

•^Deduced from the matrix element < 7(92, /3)r^''7((7i, a) >, where T^^"^ is the energy 
momentum tensor, ^ 
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1.3 The Equivalence Principle and Causality 

The equivalence principle exists in two forms: weak and strong. The weak 
equivalence principle states that at each point in spacetime there exists a local 
Minkowski frame, which is a fundamental requirement of general relativity^. 
The strong equivalence principle (SEP), on the other hand, states that the 
laws of physics are the same in all LIFs at different points in spacetime, 
and at the origin of each LIF they take the special relativistic form. Then 
the coupling of curvature to the electromagnetic field in the effective action, 
Eqn. ()1.4j) . is a violation of the SEP. Due too this violation of the SEP, QED 
in curved spacetime remains a causal theory - despite the modification to the 
physical light cone. 

This can be seen more clearly by considering Global Lorentz invariance^ 
(GLI), which is the special relativistic equivalent of the SEP. In special rela- 
tivity GLI states that faster than light signals automatically imply the pos- 
sibility of unacceptable closed-time-loops. Fig. 11.21 That is, if you can send 
a signal backwards in time in one frame, then it should be possible in any 
frame. However, if you break this GLI, a signal backwards in time in one 
frame does not automatically imply you can send a signal backwards in time 
from any frame. Therefore, in the case of QED in curved spacetime, due to 
the breakdown of the SEP we retain the fundamental property of causality: 
faster than light signals can be seen to go backwards in time in a certain 

frame, but this no longer implies that you can send a signal backwards in 

"'This implies that general relativity is formulated on a Riemannian manifold 
^Global Lorentz invariance states that the laws of physics are the same in all inertial 
frames. 
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Figure 1.2: An unacceptable closed-time-loop. 

time from any frame. Thus, in this case, spacelike motion does not necessarily 
imply a causality violation^. 

1.4 Critical Stable Orbits 

One of the simplest curved spacetimes is the (Ricci-flat) Schwarzschild space- 
time, which describes the spherically symmetric geometry outside a star. The 
geometry of such a spacetime structure is given by the line interval: 

ds-' = (1 - _ (1 _ '^y^dr^ - r^de' - (r^ sin^ e)d(j)^ (1.11) 

In this spacetime the equations of motion, derived for k"^ = 0, have a special 

solution for null geodesies. This solution describes a light ray in a circular 
^Further discussion of causality is given in QH] 
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orbit with a radius r = 3M and an associated critical impact parameter. 
Therefore, a photon projected from r — oo with the critical parameter tends 
to the stable circular orbit by spiralling around it. But, in the context of 
vacuum polarization effects, we may derive the null equation of motion us- 
ing the quantum modified Maxwell equation. In this way the orbit equation 
should give us stable circular orbits that are dependent on the photon polar- 
ization, i.e. splitting the circular orbit depending on polarization. Therefore, 
a photon coming in from r — > oo, with a particular polarization and impact 
parameter, should tend to its associated critical stable orbit by spiralling 
around it in the classic way. 



Chapter 2 

Null Dynamics in Schwarzschild 
Spacetime 

In this chapter we will give a brief overview of null dynamics in Schwarzschild 
spacetime. Starting with the geodesic equations and the line interval for pho- 
tons we will derive the orbit equation, which will then be used to determine 
the critical stable orbit for photons and the associated impact parameter. 
We will then go on to show that decreasing the impact parameter, from the 
critical value, causes the photon trajectory, given as a function of 0, to spiral 
into the singularity, but as a function of t it tends to the event horizon. The 
main aim of this chapter is to familiarise ourselves with the classical solutions 
of photon trajectories in Schwarzschild spacetime, so we can then compare 
the equivalent results for the quantum modified case. 

2.1 Equations of Motion 

When we consider motion in a plane, where ^ = | and 9 = 9 = 0, the 
geodesic equations for Schwarzschild spacetime, ()A.lflj) - (|A.13j) . can written 
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as: 



2M,.. 2M . , , 

^-—)t+—'^t = ^ (2-1) 



^2 _ ^^2 ^ ^ -1^ _ n -2-2 _ Q 2.2 



2rf0 + r^^ = (2.3) 

And the spacetime line interval becomes: 

2M, , 2M, 1 9 o ;o , , 

{1 - —)t^ - {1 - —y^r^ -r^^^ = K (2.4) 

where = for photons and ±1 for space-like or time-like motion. Now, 
to derive the orbit equation ^ we use Eqns. (j2.1j) . (|2.3jl and the interval 
equation (12.4(1 . in the simplified form: 

= t + —{l-—rHt (2.5) 
2f ■ 

= 0+— (2.6) 

r 

2M, 9 , 2M,,., 9, 2M , , 

1 = -r^+ 1 V-r^l W 2.7 

In order to specify photon trajectories we can set = at anytime, however, 
then r would be interpreted as an affine parameter and not proper time. To 
proceed we divide (jSH)) by f- and by ^: 



i 2M , 2M, -I 

T + — 1 

t r"^ r 



These can written as: 



d) 2r , , 

= -^- + — 2.9 

ch r 



= --(lnt + ln(l - — )) (2.10) 
rl ■ 2r 

= -(ln0 + -) (2.11) 
dr r 
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Then, solving these we have: 

i = (2.12) 

r 

= ^ (2.13) 

Where E and J are constants of integration denoting total energy and an- 
gular momentum about an axis normal to the plane 6 = tx 12. Now, using 
Eqns. (I2.12|l and ()2.13|) to substitute for i and in Eqn (|2.7|1 and rearranging, 
we have: 



, ^, 2M,.P , 2M, , 2M,^, 2M, , n , 
+ r\l — + K{1 = 1 2 1 2.14 



Rearranging and simplifying, this becomes: 



2 



(-)2 + (l )(i^+ ) = e2 2.15 

dr r r'^ 

,dr, , 2M,,K ,^1 , , 

(^) = J5|l-(l--)(^ + ^)l. (2.16) 

We now have the components of the four momentum P" in spherical polar 
coordinates 



= (t,r,^,0) 

= E{F^'Al-n§^ + ^)]K^.^) (2.17) 

where D = ^ is the impact parameter, P = (1 — ^) and = for photons 
and K = +1 for particles. 

2.1.1 Orbit Equation 

In order to derive orbit equations that are physically understandable we need 
to represent them as r(0), r{t) and 0(t). In this form we can analyse the 
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orbit paths as a function of rotational angle and the time taken to reach a 
certain point along the angle. Therefore, by using 



dr dr dcj) dr J 
dr dd) dr dchr"^ 



(2.18) 



we can rewrite (I2.15|) as 



dr P 9M P 

= Qr- + ii-—)iK+'-) 

dq) r 
^ (_)2 = ^E''-K)— + 2MK—-r^ + 2Mr (2.19) 

Now, transforming u = r^^, so at r = cxd we have u = 0, Eqn. ()2.19|1 becomes: 

(^)V^ = ^E'-K)^ + 2MK^-r' + 2Mr 
do 



,dU,r, (E'^—K) M O , , Q 

{-^Y = - — j2 — ^- + 2MK—-u^ + 2Mu^ (2.20) 



Doing similar manipulation for (pit) we have 

d(j) d(j) dr 
dt dr dt 



(2.21) 



and using Eqns. p.l2|) and ()2.13j) in this, and transforming m = we have: 

dd) D, 2M, o, . , . 

-£ = -{1 - —) = Du\l - 2Mu) (2.22) 

Finally, we can determine u{t) by inserting: 

du du dt du,^ ^, , ^ X H 1 

— = - = —[Du\l - 2Mu)r^ 2.23 

dd dt d(h dt^ ^ ^ ' 



into Eqn. ()2.20j) . which gives 



[^f = [Du\l - 2Mu)f[ ^^\J^^ + 2MK^ - + 2Mu^] (2.24) 
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Now, the Eqns. (|22DI), (l2T22ll and (12:2411 are the general equations that de- 
termine photon and particle trajectories in the plane B = t\ jl. If we take 
= we then have the required photon trajectory equations 

= -^-«^ + 2M«3 = /H (2.25) 



{Du^{\-2Mu)ff{u) (2.26) 



[^f = [Du\l-2Mu)]^[—-u^ + 2Mu^ 



We can also write another equation, specifically for null radial geodesies. By 
using the fact that = 0, Eqn. (j2.13l) implies J = 0, then Eqns ()2.12jl and 
(|2.15j) become 

2M,dt 

1 V = ^ 2.27 

r dr 

(IT 

^ = ±E (2.28) 
dr 



Combining these we have: 



dr , 2M, 

^ = ±(1--) (2.29) 



2.2 Orbits in Schwarzschild Spacetime 



We will now discuss the solutions of Eqns. (I2:25ll . (12:261) and 1^2^. We will 
start off with the simple case of the radial geodesic and then go onto the case 
of the general orbits. For the general case we will solve ()2.25|] and ()2.26j) 
to determine the critical stable orbits and the associated impact parameters. 
We will then go on to study the trajectories of photons as they fall into the 
singularity 
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Figure 2.1: Radial Null geodesies take an infinite coordinate time to reach 
the event horizon at 2M 



2.2.1 Null Radial Geodesic 

Eqn. ()2.29|) can be solved by rewriting it as: 
dt , 1 

- = ±1 2.30 

ar r 

=^ t = ±{r + 2M \og{^ - 1)) + constant± (2.31) 

This solution (specifically the — one) gives the trajectory of a photon coming 
from infinity and into the black hole. It shows that the photon takes an 
infinite coordinate time to reach the horizon, which can be seen in Fig. 12.11 
However, if we solve Eqn. ()2.28|) . i.e. in terms of the affine parameter, we 
can show that the photon reaches and crosses the event horizon without ever 
noticing it, 

r = ±Et + constant± (2.32) 
this can be seen in Fig. 12.21 
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T ZM 




r 



Figure 2.2: Radial Null geodesies reach and cross the event horizon at 2M 
without noticing it when defined with the affine parameter 



2.2.2 General Null Geodesies and Critical Orbits 



We will first solve Eqn. ()2.25|) in order to determine the the critical stable 
photon orbits. In order to do this we first consider the point of equilibrium 
and the associated impact parameter. This equilibrium point occurs when 

du 

which means as the photon is orbiting the black hole u does not change; so 
if the radial distance does not change it implies a circular orbit. Therefore 
we must solve 

j{u) = l^-u^ + 2Mu' = (2.34) 
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The sum and product of the roots Ui, U2, and M3 of this equation are given 

Ui + U2 + U3 = — and ^1^2% = - 2MD^ 
This shows that f{u) = must have a real negative root, and the two 
remaining roots can be real (distinct or coincident) or be a complex conjugate 
pair; however, the occurrence of coincident positive real roots implies the 
existence of a circular orbit. Thus, if a coincident root occurs it should be at 
the point given by the derivative of f{u) 

f{u) = 6Mu'^ -2u = (2.36) 

Which then has the solution ui = U2 = (3M)~^. For this solution the impact 
parameter of equation ()2.34|) is -D = (3\/3)M. From the product condition, 
equation ()2.35|) . we find that the roots of /(n) = are 

Ui=U2 = ^ and % = --^ and D = {3V3)M (2.37) 
6M blVI 

Therefore, when the impact parameter is -D = (3v^)M then ^ vanishes for 
u = (3M)~^, which implies a circular orbit of radius 3M is an allowed null 
geodesic [2]. 

Now we can consider a photon at u = with an impact parameter D = 
(3v^)M. This, then, gives a trajectory of a photon spiralling in and tending 
to the critical orbit at n = (3A/)~^. The general differential equation for 
this impact parameter is given by rearranging and substituting for D in 
Eqn. 

(|)^^2M(„,.gL)(.„__L)3 (2.38) 



^for au^ + bu^ + c = we have ui + tt2 + U3 = — ^ and U1M2M3 



CHAPTER 2. NULL DYNAMICS IN SCHWARZSCHILD SPACETIME 18 
From [2| we have the solution to this as 

where 0o is a constant of integration, given by: 

tanh2(-i0o) = ^, (2.40) 

which gives: n = (r ^ oo) when = 0, and u = when (f) oo. 
Therefore a null geodesic arriving from infinity with an impact parameter 
D = (3v^)M approaches the circle of radius 3M, asymptotically, by spi- 
ralling around it, as can be seen in Fig. \2.3f . Also, numerically solving 
Eqn ()2.26|) we can show that as time increases u tends to which can be 
seen in Fig. 12.41 

Finally we can show that when an impact parameter other than D = 
(3-\/3)M is used, for example if we set D = (3-\/3)M — 0.1, the solution 
of Eqn. (j2■25^) shows that the photon falls past the critical orbit, through 
the event horizon, and into the singularity. Fig. 12.51 Using this new impact 
parameter in Eqn. ()2.26|) and, again solving numerically, we see that the 
photon comes in from infinity and tends to the event horizon, u = 
asymptotically in time t, Fig. 12.61 



^This figure was plotted using Eqn Ij2.39ll . We also obtained the same plot by numeri- 
cally solving Eqn H2.38II in Mathematica. 
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u = l/3M 



— — ^ — . — . — . — . — . — . — . — . — . — . — . — . — . — . — . — . — . — . — . — - 9 

£ 4 f $ 10 

Figure 2.3: Null geodesic, with impact parameter D = {3^/3)M, arriving 
from infinity and approaching m = 3^ asymptotically (M=l/3) 




U=1/3M 



£0 40 £0 SO 100 120 140 

Figure 2.4: For impact parameter D = (3a/3)M a null trajectory tends to 
u = ^ asymptotically with time. 
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Figure 2.6: Null geodesic, with impact perimeter D = (3-\/3)M — 0.1, arrives 
from infinity and tends to m = ^ asymptotically in time t. 



Chapter 3 

Quantum Gravitational Optics 



3.1 Photon Propagation in Curved Spacetime 

Maxwell's equations, in curved spacetime, 

= D^F^"" (3.1) 

= D^F,x + D,Fx^ + DxF^, (3.2) 
F^, = d^A, - d,A^, (3.3) 

cannot be solved explicitly, and even in cases of extreme symmetry explicit 
solutions are difficult; this is due to the fact that as curved space acts as 
a dispersive material (i.e. bending light rays) plane wave solutions do not 
exist. 

3.1.1 Geometric Optics in Curved Spacetime 

In dispersive materials, where light rays are bent, we can consider the solution 
of Maxwell's equations to be a simple perturbation of the plane wave solution. 
For example in curved space, relative to an observer, the electromagnetic 
waves can appear to be plane and monochromatic on a scale that is much 

21 
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larger compared to the typical wavelength, but very small compared with 
the typical radius of curvature of space time. Such "locally plane" waves can 
be represented, in geometric optics, by approximate solutions of Maxwell's 
equations of the formp]: 

J^^" = Re{Fr + teFf + . . (3.4) 

and the electromagnetic field vector, defined by Eqn. (I3.3jl . takes the form: 

= Re{A>; + isA''^ + . . .)e'' (3.5) 

where the electromagnetic field is written as a slowly- varying amplitude and a 
rapidly-varying phase. The parameter e is introduced in order to keep track 
of the relative order of magnitude of terms, so in curved space Maxwell's 
equations can be solved order-by-order in e. In this formulation the wave 
vector is then defined as the gradient of the phase of the field, = Dfj_^9 = 
d^^9, which in terms of the quantum interpretation is identified as the photon 
momentum. We can also write A'^ = Aa'^, where A represents the amplitude 
and a^, (normalized clS Qj^Qj — —1) specifies the wave polarization. These 
vectors then satisfy the condition k^a^ = 0. 

Geometric Optics and Null Dynamics 

In this notation Eqn. ()3.H) can be written, to leading order 0{^), as: 

^ k^Fr = (3.6) 
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and Eqn. ()3.3|1 becomes: 

= '-{d^eA^,-d,eA^,)e''- 

F^^^ = k^Ai^-KAi^ (3.7) 

Now, by combining these, we have: 

k,Fr = - = 

= (kfj_k^)A'^ — k'^(k^A^) 
^ k^a" = (3.8) 

and from this we can deduce that k"^ = 0, i.e. k'^ is a null vector. Also, it 
follows from the definition of k'^ as a gradient that D^k,y = D^k^, so 

PD^r = k^D^kf, = ^D'^k^ = (3.9) 

Using this, and the fact that light rays are defined as the curves given by 



X I, 



s) where ^ = /c^, we can derive the geodesic equation as follows [Z|: 



= k^'D^k'' 

d'^x" ^„ dxf^ dx^ 

3.2 Quantum Modified Null Dynamics 

As was seen in Sec. II. 2[ using the effective action, Eqn. (|1.4p . the equation 
of motion becomes: 

= - ^[2hR^^D^^F'^- + Acg^^R^^^pD^'F'^P] (3.11) 

and the Bianchi identity, Eqn. ()3.2|) . remains unchanged. Now, as before we 
can determine the quantum modified light cone and geodesic equations. 
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Quantum Modified Light Cone 

Substituting Eqn. dSH in (ITTT^ we find, again to 



= -iD,e)e^lFr - -^[2bR,,-iD'^e)e^l Fl 
e m'^ e 



+ Acg-^R,,^p'-{D^^e)e'TF^P] 

^ = A;^Fr - \ [2hR^^k^F^'' + Acg"^ R^^^^k^ F^'] (3.12) 

and using Eqn. (|3.7jl we have: 

{) = k^{k^^A'i-k^A'i) - \[2hR^y^k^^{k^A'i-k''A'l) 

+ Acg'^R^rXpk^ik^Al-kPA^,)] 

using k^A^ = this becomes: 

= k^k''A'( - l^[2bR^xk^ik'^Al-k''A^,) 

+ 4cg''^R,rXpk^{2k^Al)] (3.13) 

where A'^ = Aia^ , and the last line is simplified by relabeling of indices. 
Now, Contracting with Aa,y and eliminating A we have: 

= ki^k^ H -Rn\k^k -RnrXpk^k aiCii 

This then gives the quantum modified light cone: 

- ^R^k^k^ - ^R.rXpk^k^a^aP = e - 5k(a) = (3.14) 

where we have replaced the constant c = with Cq, = 3^ for convenience 

of interpretation, i.e. the sign of the light cone is immediately obvious from 
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R^rXpk^k^oT' . As we will be working in the Schwarzscliild spacetime the 
quantum modified light cone for the Ricci flat case {R^\ = 0) is given by: 



e = -^R^^^pk^'k^a^aP = 6k{a) (3.15) 



Here the sign of the light cone depends on polarization and the photon tra- 
jectory; if the correction is positive we have space-like motion, and if it's 
negative we have time-like motion. 

Quantum Modified Geodesic Equation 

The photon trajectories corresponding to the quantum modified equation of 
motion, (jH.llj) . can be represented by a generalised version of Eqn. (j3.9j) : 



= ^D^e - ^R.xk'-k^ -^R.^.pk'-k^a^a'' 
2 mi 



= -D,k^ - —D^bRf^xk^k^ + 4caR^rXpk^k^a^a''] 



2 ^ 

<fx^ dx^ dx^ 



= ^^ + r 



ds^ ds ds 

- —d^libRfs^ + ic^Rp^^^a'^a^) — —] (3.16) 

where we have used k^^ = and covariant derivative in the second term is 
replaced by a partial derivative as it's acting on a scalar. In Ricci spacetime 
this equation becomes: 

3.3 Horizon Theorem and Polarization Rule 



There are two general features associated with quantum modified photon 
propagation[2j. First, it is a general result that the velocity of radially di- 
rected photons remains equal to c at the event horizon. Second, for Ricci 
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flat spacetimes (such as Schwarzschild[T| and Kerr|H] spacetimes), the veloc- 
ity shifts for the two transverse polarizations are always equal and opposite. 
However, this is no longer true for non-Ricci flat cases (such as Robertson- 
Walker spacetimejl]). In these cases, the polarization averaged velocity shift 
is proportional to the matter energy-momentum tensor. These features can 
be easily shown by using the Newman-Penrose formalism: this characterises 
spacetimes using a set of complex scalars, which are found by contracting 
the Weyl tensor with elements of a null tetrad[2]. 

Newman-Penrose Formalism 

We choose the basis vectors of the null tetrad as[Z|: = k^, the photon mo- 
mentum. Then, we denote the two spacelike, normalized, transverse polar- 
ization vectors by and 6^ and construct the null vectors = -^{a^ + ib^) 
and rh^ = "^(^^ — ib^)- We complete the tetrad with a further null vector 
n^, which is orthogonal to m'^ and m^. We then have the conditions: 

l-m = l- rh = n- m = n- m = (3.18) 
from orthogonality, and: 

l-l = n- n = rn-m = rh-rh = (3.19) 
since the basis vectors are null. Finally, we impose: 

I ■ n = —m ■ rh = 1 (3.20) 
The Weyl tensor, given in terms of the Riemann and Ricci tensors, is: 



1 



C^i/'yS R^vyS Q i.Vfi'f'^i^s ^w^fR^s V^J.S'l^u"/ ~l~ ^usR^ 



V-1) 



2 

QiVf^-rVi^s - 'n^i5r]u-^)R (3.21) 
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where i?^^ = rj'^^R^u'yS and R = rj^'^R^y] and the Weyl tensor satisfies the 
trace-free condiation: 

V^'C^y^s = (3.22) 

and the cyclicity property: 

Cl234 + C'l342 + C'l423 = (3.23) 

Now, using the null tetrad, we can denote the ten independent components 
of the Weyl tensor by the five complex Newman-Penrose scalars: 



^0 = 




^1 = 




^2 = 




^3 = 




\i/4 = 





(3.24) 

3.3.1 Polarization Sum Rule 
Ricci Flat Spacetime 

In Ricci flat spacetime, summing the quantum correction over the two po- 
larizations leads to the following polarization sum rule: 

^5A;(a) = (3.25) 

a 

This can be proven by suming the quantum correction in Eqn. ()3.15|1 over 
the two polarizations, 

5k{a) = ^Y^ Rf^u^sk^k-^a'a^ (3.26) 
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In the Newman-Penrose basis, using k = I, a = ^{m+m)^ h = —'^{m—rh), 
and the fact that Cf^^-yS = Rfiu-yS for the Ricci flat case, we have: 

a 

= C^^^sl^Pim^rh^ + m''m^) (3.27) 

This particular contraction is equal to zero as it's not part of the complex 
scalars in Eqns. ()3.24|) : hence the sum of the two quantum corrections is zero. 
This implies the trajectory (and velocity) shifts are equal and opposite. 

Non-Ricci Flat Spacetime 

For the non-Ricci flat spacetimes the polarization sum rule is given as: 

Sk(a) = -^(2b - 8c^)T„^Fr (3.28) 
mi 

a ^ 

where T^^, is the energy-momentum tensor. This can be shown by proceeding 
as before, but now we include the Ricci tensor and scalar, as in Eqn. (|3.2H) : 

^ R^^^sk^'k^d'a^ = C^^^sl^Pim^m^ + rffm^) - R^^FP (3.29) 

a 

As before the first term on the RHS is zero, and the second term is only 
dependent on the photon momentum. Then, combining this with the Ricci 
term in Eqn. ()3.14|) we have: 

6k{a) = i-(26 - 8ca)R.,k''k' (3.30) 

a ^ 

Finally replacing the Ricci tensor with the energy-momentum tensor, by 
using the Einstein equation 

Rf,u = -SnT^u + ^Rg^u (3.31) 
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we obtain Eqn. (ITM . 

3.3.2 Horizon Theorem 

At the event horizon, photons with momentum directed normal to the hori- 
zon have velocity equal to c, i.e. the light cone remains k"^ = 0, independent 
of polarization^. 

This can be easily proven for the Ricci flat spacetime, using the orthonormal 
vectors k" = {Ek,Ek, 0,0), a'' = (0,0,1,0) and a'' = (0,0,0,1). Therefore, 
using these vectors in Eqn. (|3.15|) . we have: 

e = ^Rabcdk^k^a'^a'^ = (3.32) 

So, in Ricci flat spacetime all radially projected photon trajectories remain 
unchanged. 

It is also possible to prove the horizon theorem for the general case (for 
Ricci and non-Ricci flat spacetimes) that the light cone at the event horizon 
is unchanged. This can be seen in the null tetrad, so that the physical, space- 
like, polarization vectors and lie parallel to the event horizon 2-surface, 
while fc^ is the null vector normal to the surface. Then, from Eqn. ()3.14j) we 
have for the two polarizations: 

k = ~Rn'yk^k^ ~t~ — Ri^jy^gk^ k^ CL a 




Using Eqn. (I3.31|) and the fact that Cf^u-yS^P {fn'^fn^ + m'-'m^) = 0, we can 
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write this as: 

^2 ^ _^^2b - Ac^)T.J''P ±^C^,^sl''Pirn''m^ + rh'rh^) (3.34) 

and in terms of the Newman-Penrose scalars this can be written as: 

fc2 ^ _^(26_4c„)r^^Pr±^v|/o, (3.35) 

where the simphfication in the last term on the RHS is possible as \E'o is 
real for Schwarzschild spacetime. In general Eqn. ()3.3■5^) is non zero, how- 
ever, at the event horizon the terms: T^^l^l^ and \E'o are zero for stationary 
spacetimes[3 ini^- 

Physically the Ricci term represents the flow of matter across the horizon 
and the Weyl term represents the flow of gravitational radiation[2j, and both 
are zero in classic general relativity; and as, even with the quantum modi- 
flcation, the light cone remains unchanged at the event horizon, means the 
event horizon is fixed and light cannot escape from inside the black hole. 



^Stationary spacetimes are independent of time and may or may not be symmetric 
under time reversal 



Chapter 4 

Quantum Modified Trajectories 



In this chapter we will analyse how the classical null trajectories, described in 
Chapter El are modified in Schwarzschild spacetime, due to quantum modifi- 
cations of the equations of motion of general relativity. Using Eqn. (j3.14j] we 
will calculate the quantum corrected version of Eqn. (j2.25j] . which will then 
describe the quantum modified motion of a null trajectory in a Schwarzschild 
spacetime. Using this, and following the methods of Chapter O we will begin 
by studying simple critical circular orbits at the stationary point, u = 1/3M. 
As stated in the polarization rule, the critical orbit should be shifted up and 
down by equal amounts, depending on the polarization of each photon. Also, 
we will show that these modifications are only valid if the "classic" impact 
parameter corresponding to the stationary orbit is adjusted, depending on 
the photon's polarization, in order to compensate for the orbit shift. This 
information, of the modified orbits and the corresponding impact parameter 
for each polarization, will then be used to determine the general trajectory 
of a (vertically or horizontally polarized) photon coming in from infinity and 
tending to one of the two shifted critical orbits. We will then go on to show 
that these quantum modifications have no effect on the event horizon, that 

31 
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is, when the impact parameter is accordingly adjusted and a photon falls into 
the singularity the event horizon remains fixed at m = 1/2M. 



4.1 Quantum Modified Circular Orbits 

In Schwarzschild spacetime, using k'^ = dx^'/dr, Eqn. ()3.15|) can be written 
as: 

= - (1 - — )V + r2(l - — 

+ {l-—y^' + ^{l-—)R,.,skWa' (4.1) 



and Eqn. (ITTTI) as: 



.. 2M, 2M, 1 . 
= t + —{l-—y^rt 

- ^m.^.sa'^a'^^) (4.2) 
mj dr dr 



r mf dr dr 



= ^ + — -^d^{R^,-ysa''a' ——) 4.4 
r mf dr dr 



= -^MR^.u^sa'a'^^) (4.5) 



e 



where = Now, in order to consider quantum modified circular orbits 
we require three things: (i) the Riemann tensor components, (ii) the photon 
wave vectors, k^^ = for circular orbits, and finally (iii) the photon polar- 
ization vectors, a'^. Due to the circular nature of the orbit the simplest basis 
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to work in is the orthonormal basis. In this frame the required polarization 
and wave vectors, for circular orbits, are simply given as: 

= (0, 1, 0, 0) Planar Polarized (4.6) 



a'' = (0, 0, 1, 0) Vertically Polarized (4.7) 

= —— = {Ek, 0, 0, Ek) 4-momentum along (4.8) 
ar 

Now using these photon vectors, the quantum modifications in Eqns. 1)4.111 - 
()4.oj) . can be expanded to give: 

• For the planar polarized case we have: 

Rabcdk k a a = Rarcrk k = Ef^i^R^rcrk k ) 

= El{Rtrtr + R(l>rcf>r + '^Rtripr) (4-9) 

• For the vertically polarized case we have: 

Rabcdk^k^a'a" = R^ecek'^k'' = El{Raecek''t) 

= El{Rtete + R(t>e<f>e) (4-10) 

Using the six independent components of the Riemann tensor in the orthonor- 
mal basis^, we find: 

Rtrrt — —Rtrtr — =^ Rtrtr — (^-H) 



^ Which we have calculated in Appendix A 
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From Eqns. ()3.15|) . ()4.16|) and ()4.18|1 we can see that: 



k'^ ~ — \ — For Planar Polarization (4.20) 

,o 'iElM 

~ \ — For Vertical Polarization (4.21) 



This implies that, as the light cone for planar polarization is positive, it 
represents a photon trajectory with a speed < c, and as the light cone for 
vertical polarization is negative, it represents a photon trajectory with a 
speed > c. 

Planar Polarization 

Considering the planar polarized case first, we can use Eqns. ()4.16j) and ()4.17p 
to rewrite Eqns. as: 

= fe + r^{i 

2M, 8c^, 2M^3ElM 
+ 1 + 4.22 

= t + —{l-—)-'rt (4.23) 



M , 2M, 1 . , 2M, -p 
= f-— (1-— )^V2-r(l 



M, 2M,.. 4:Cc,9ElM 
+ — 1 )t' + ^^- 4.24 



2r0 

+ — (4.25) 



where Eqn. ()4.5jl becomes zero. Now, using the solutions of (|4.23j) and (jliini), 
as given by (|2.12j) and (|2.13|1 . we can rewrite Eqn. ()4.22|1 as a simple quantum 
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modified trajectory equation for circular orbits with radius r and in a plane 



9 = 

u 2 • 



r-E^ + (l )— H ^(1 ^ " 



^9 ,dr^^ , 2M^,2AcElM J\ 
- = + (4.26) 

To make this equation more meaningful and easier to solve we make the 
following transformations^: 



dr d(j) dr J 
d(f) dr d(f) 



1 du dr du , 1 ^ , . ^ 



Therefore, Eqn. ()4.26|1 becomes: 



Which can be written as: 



(4.27) 



■ (- A)' ■ {Ju'? + (1 - 2M«)(^^^^^4^ + '^'«') 

(XCp U TDj^ 

)2. J2 + (i_2Mm)( ^4 + JV) 4.29 

dm mi 



'^f = 2M^Au^ + {2M-MAy-u^ + ^ (4.30) 



where we have defined the dimensionless constant: 



(4.31) 



^We must note that i? and Ek are not equal, is the energy from the classic relativistic 
orbit equations, while Ek is the quantum energy of the photon 

dr 



^Using J = r2^ 
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In the last form we have used (without proof) the relation = a/S-E, which 
will be proven in Sec. 14.2.1 L Eqn. (|T68ll .'^ 

Vertical Polarization 

Similarly, for the vertically polarized case, using Eqns. ()4.18j) and ()4.19|1 in 
Eqns. (PTT|l - (p3|l to give: 

= - (1 _ _)V + r2(l - — 

+ (1 )rV-^(l (4.32) 



t + — (1-— )-Vt (4.33) 



= r-— 1 rV^-rl 02 

+ — 1 ^ -^^^ 4.34 



= + — (4.35) 

r 

which, in a similar way as before, gives us the equation of motion for the 
vertically polarized photon: 

rfw 2 ^ -2MMm^ + (2M + MA)u^ - + ^ (4.36) 

4.1.1 Quantum Modified Critical Orbits 

Now the general equation for the quantum modified circular orbits is: 

A)2 ^ ±AM(2Mu - l)u^ + 2Mu^ - + 2^ = f(u) (4.37) 
ad) -D 



'^Ca is dimensionless, while D and me have dimensions of length and inverse-length 
respectively 
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where + is for planar polarization in the r direction, — is for vertical polariza- 
tion in the 9 direction and D is the impact parameter. As A — ^ Eqn. ()4.37p 
tends to the classic orbit equation in general relativity, Eqns. ()2.25j] . There- 
fore, in order to determine the magnitude of the quantum correction we can 
calculate the order of A using typical values for mg, D and Cq, in Eqn. ()4.3ip . 
Using^ niec/h ~ 10^^ for the electron mass (as it is given as inverse length), 
Ca = a/SGOir ~ 10^^, and the mass of the sun inserted into the critical im- 
pact parameter: D = 3 a/3G Mq/c^ ~ 10^ (given in terms of length), this 
then gives us^: 

^ _UCa iU 32 . . 

m2D2 (1011102)2 

With the order of A being so small the correction in Eqn. ()4.37|) will be tiny 
compared to the size of the orbit (r = 3GM/c^ ~ 10^); therefore the modified 
orbits will not differ from the classic critical orbit, given by Eqn. ()2.25ll . by 
very much. 

We will now determine the quantum modified critical circular orbits. In 
order to solve Eqn. (|4.37jl . we can use the fact, from the polarization rule, 
that the modified orbits should be shifted above and below u = hj equal 
amounts depending on polarization. So we expect a solution of the type 
u = gj^i^M; which means we can try a simple modified solution of the form: 

u = Uq + kui (4.39) 

where Uq = g^j, ui is the quantum modification, and is a small constant 

^As we were working with G = c = h ^ 1, we must reintroduce these constants to 
obtain the correct order of A 

^This result is also shown in 
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depending on the quantum modification A. 
Planar Polarized Critical Orbit 



Working with Planar polarization, we can substitute the solution ()4.39j) into 
the derivative df /du of Eqn. ()4.37|) : 

^ = 2AM\^ + 3AM(2Mu - l)u^ + QMu^ - 2u = (4.40) 
du 

and as A <^ 1 and k A, then the only terms of relevance are the ones first 
order in k and A, everything else can be assumed to be ~ 0. Therefore, we 
have^: 

2AM^ul + QAM^ul - 3AMul + QMul + 12MkuiUo - 2uo - 2kui = 



AM , 

^ ^1 = -^«o (4-41) 
Now, substituting this into the trial solution (I4.39jl . we have: 

AM , 

u = uo{l + —ruo) 4.42 

D 

where mq = g^j is the classic solution. Therefore the classic orbit uq is 
modified by to first order in A. Also, with this orbit modification we 
require an associated, modified, impact parameter, which should take the 
form: 

1 1 

where 1/-Dq = 1/27M^. The modified impact parameter can be found by 
substituting (|4.42|] and (|4.43jl into (I4.37jl and solving for /?. Doing so, we 



'u — uq + kui u — Uq + 2kuiUo u — Uq + Skuiu^ 
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find 



8. 



[^f = 2M^Au'-AMu^ + 2Mu''-u' + ^ + f] = 
ad) Dn 



2M^Aul - AMul + 2Mul + MMm; 



""o 



- 3^^^^o + ;^ + /? (4.44) 

We have 2Muq — Uq + = 0, as this forms the classic equation of motion 
for circular orbits. Therefore, 

2M'^Aul - AMul + M'^Aul - -AMul + P = ^ 

3 



Now, as Mo = -TTf, we have: 



(3 = ^ul (4.45) 

Therefore the modified impact parameter, for planar polarization, is: 

1 1 AM ^ 

ln = Wi^ — "» 

Then, substituting for Dq, we have: 

1 _ 1 AM 3 _ ul AM 3 
~ 3(3M)2 + ^""o ~ y ^ ^""^ 

= ^(l + AMuo) (4.47) 
Vertical Polarized Critical Orbit 

Doing the same for the vertically polarized photon, i.e. by using: 

df 

-i- = -2AM'^u^ - 3AM(2Mu - l)u'^ + QMu'^ - 2u = (4.48) 
du 



^we, again, work to first order in A: u = uo + A^Uq, = Uq + A^-u^, — + A^-Uq 
and u*^ul + A^ul 
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and substituting the trial solution ()4.39|) we find: 

-2AM^ul - 6AM\l + SAMul + 6Mul + UMkuiUo - 2uo - 2kui = 



AM , 

^1 = --^4 (4-49) 



AM , 

u = no(l —uo) (4.50) 

which is equal, but opposite in sign, to ()4.4H) . as is expected from the po- 
larization rule. Also, as before, the relevant impact parameter is given by 
substituting ()4.5fl|) and ()4.43|) into the negative equation of ()4.37|) and work- 
ing to first order in A. 

^ = -2MMu^ + AMu^ + 2Mu^ -u'^ + ^ + (3 = 
d(p Dq 

= -2M^Auq + AMul + 2Mul - MAu^ - ul 

1 , 1 

\ 

Eliminating terms and rearranging, as before, we find: 



^AMul + —^+^ (4.51) 



AM n 

(3 = -^ul (4.52) 

Therefore the modified impact perimeter, for vertical polarization, is: 

1 1 AM ul AM , 



2 Q/'QA/f^2 Q "0 



L)2 3(3M) 



Un = — — Mn 



?y2 

^(1-AMwo) (4.53) 



The Modified Orbits 



We now have the circular orbit solutions for Eqn. ()4.37|) and the relevant 
impact parameters: 
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Figure 4.1: Orbit shifts about the classic critical orbit Uq = 1/3M, with 
corresponding shifts in the speed of light (Not to scale) 

• Planar (r) polarization solution (c<l) 

M \ 

u = u^^A{—ul) — = ^(1 + AMmo) (4.54) 

• Vertical(6') polarization solution (c>l) 

M 1 

u = u^-A{—ul) — = -^{\-AMu^) (4.55) 

which are displayed in Fig. 14.11 We can note that as the constant A is of the 
order 10"'^^ these modifications are extremely small. 

4.2 Quantum Modified General Geodesies 

Using the solutions for the critical circular orbits and the associated impact 
parameters we can now study how the general photon trajectories are modi- 
fied due to quantum corrections. In the classic case when a photon comes in 
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from infinity, with the critical impact parameter, it tends to the critical orbit, 
as was shown in Fig. 12 .3^ and if we slightly decrease the impact parameter 
the photon spirals into the singularity. Fig. 12.51 We will now construct a gen- 
eral quantum modified equation of motion and determine how the geodesies 
change, depending on polarization, as they tend to the critical orbits. 

4.2.1 General Vectors 

From the quantum modification term in Eqn. 1)4.11) we can see that in order to 
construct a general quantum modified geodesic equation we require general 
photon polarization and wave vectors. The wave vector, k^, can no longer 
be represented by a simple constant vector pointing in the direction, as 
given in the orthonormal frame. And, even though the vertical polarization 
will remain a constant, as before, the planar polarization will now be a more 
general vector, constantly changing as the photon moves through the plane. 

General Wave Vectors 

Our general wave vector will be of the form: 

(4.56, 

which respects the light cone condition: 

= g^^k^k^ = FP - F~^r^ - r^'^ (4.57) 

where gf^u is the Schwarzschild metric. Using previous results, of Eqns (I2.12j) 
and ()2.13j) . and the fact that we are working in the plane, ^ = |, we can 
represent three of the wave vector components as: 
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where we have defined: 



:i - ^) (4.50) 



Now using Eqn. ()4■57^] we can write the final component as: 



where we have used D'^ = Now the general wave vector can be written 



as: 



Also: 



F = (EF-\E(l-:^)i,0,4) (4.61) 



k, = {E,-F-'Eil-^)K0,~J) (4.62) 



We now have: 



k,k^ = E^F^^ - F^^E\l - ^) - ^ = 

as required. Therefore, eliminating E^ we have the general wave vectors: 

= ^(^-^ (1-^)^0,^) (4.63) 
K = ^i^.-F~\l-^)K0,-D) (4.64) 

We now need to normalize these vectors so, as they come in from infinity 
and tend to the critical circular orbit, they correspond to the critical orbit 
vector ()4.8|) . However, as ()4.8|) is given in the orthonormal basis, and we 
are now working in the coordinate frame, we must use the tetrad ()A.44jl . 
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given in Appendix B, to transform ()4.8j] to its coordinate frame equivalent. 
Therefore, ()4.8|) in the coordinate basis is given as: 








E, 






V -r I 



(4.65) 



This now corresponds to a wave vector for a circular orbit in the coordinate 
frame. To represent the critical orbit we simply substitute for r = 3M, in 
which case F = 1 — ^ — ^ i, thus ()4.65|) becomes: 

^ = ^,(73,0,0,^) (4.66) 

Now, if we evaluate ()4.64jl with r = 3M and D = 3^/3M we have: 

F = ^v^(v^,0,0,^) (4.67) 

Therefore, ()4.67|) is similar to ()4.66j] up to a constant of normalization, given 



as: 



Ek = ^Vs = eVs 

Then, the normalized form of (|4.64|) is given by using ()4.68p : 

, ,, Eh , „ 1 , D^F .1 D . 
p = ^ F-i 1 ^ ^0,- 



General Polarization Vectors 



(4.68) 



(4.69) 



Now we need to construct the, planar and vertical, polarization vectors, a'^, 
of the photon; which must be spacelike normalized as: 



-1 a^F = 



(4.70) 
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As before, as we are working in a plane, the vertical polarization, a^, will be 
a constant, and can simply be written as: 

< = (0,0, 1,0) (4.71) 

To normalize this we do as follows: 

ai^ = g^Xi = -r\0,0,l,0) a^a^^ = -r^ (4.72) 

Therefore, the normalized vertical polarization vector is given as: 

< = (0,0,^,0) 

ai;. = -r2(0,0,-,0) (4.73) 

r 

These now satisfy both the conditions in ()4.70p . The planar polarized vector, 
flji is given in the plane of r and 0: 

a^ = (0,AO,S) (4.74) 

Now, using the two conditions in ()4.70|) we can determine A and B. From: 

= ^M.«2 = (0, -AF-\ 0, -r^B) (4.75) 

and the vector ()4.69|) we have: 

k^a2, = (1 - ^)HaF~') + {^)ir'B) = (4.76) 

a^^aa;, = A^F-^ + B\^ = -1 (4.77) 
Solving these for A and B: 

F-\ D'^F^i , 
i? = ---(l--— )M (4.78) 
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Therefore, we have: 



A = (4-79) 

(1 + ^(1-^)) ^ ' 



. DF 1 D'^F.i 
A = B = —{1 —2 4.80 



and the planar polarization vector becomes: 



a^ = (0,^,0,--(l-^)^) (4.81) 



We now have the required polarization vectors: 

1 



^(0,0,1,0) (4.82) 

ai^ = -r2(0,0,-,0) (4.83) 
r 

= ho,DF,0,-{l-^f^) (4.84) 

«2M = (0,-^,0,(1-:^)^) (4.85) 



where subscript 1 and 2 are vertical and planar polarizations respectively. 
These now satisfy the conditions in (I4.7()jl with the wave vector (j4.69j] . 
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4.2.2 Quantum Modification 

Having derived the general polarization and wave vectors, we can see from 
Eqn. ()3.15|) that the quantum modification given by: 



Sk(a) = -^Rabcdk''k'a''a'^ (4.86) 

also requires the Riemann tensor components in the coordinate frame. In 
Appendix A we have calculated the required components as: 

Using this information we will now determine the form of the general quan- 
tum modification; and from the polarization rule, this quantum correction 
should satisfy the condition: Skf = —dk^ for the two polarizations. 

For vertical (6*) polarization we have, by using (|4.69|) and ()4.82j) in (|4.86p : 

Rabcdk"'k'^o!lO'i = Radcd^k'^k'^ = Rtm^k^k^ + RrBrd^k^k^ + Ripe^d^k'^k'^ 



MF I F-^E, F-^Ek , 



MF-^.l Jl-^E, Jl-^E, 
A DEk , , DEk , 
MElF-^ MElF~\l - ^) 2ME|D2 



Sr^ 3^3 3j-5 

MElD"^ 



(4. 



Similarly, for planar (r — plane) polarization we have, by using (|4.69^) and 
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m 

~\~ k^ (x^(x<2 ~\~ fl^f^ffpk k Oj^clq^ ~\~ R.j-(j)'p(jyk k cl^Qj^ 

2M F-'Ek , DF , MF-^ DEk ^ DF ^ 



n2p 

MF-\DEk,,Ek L D^F.,DF. " 



MF .F-'Ek,2,\/^~—,2 



r y/3 



+ 



r V3 



MElD"^ 



(4.89) 



Therefore, the quantum modifications, k"^ = 6k{a), for the two polarizations 
ai and 02 (vertical and planar respectively) are: 

Sk{a,) = -(^) ^^P' Vertical (4.90) 
Sk{a2) = (%) ^^P' Planar (4.91) 

where — 5A;(ai) = 5A;(a2), as required. Now, using fc^A;^ — Sk{a)=0 we can 
write the general equations of motion for the two polarizations: 
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As before, substituting for t, and transforming r — > ^, we have: 

(55)==Bi-^"^±(— (4.92) 
We can also write m as a function of time: 

(f )^ = {Du^F)\^^ - Fu- ± {^^)Fu') (4.93) 

where we have used the substitution for A given in equation (14.3111 . Now, 
the Eqns. 14.921 and 14.931 are the general equations of motion, + for vertical 
(6*) polarization and — for planar (r-0) polarization. In these equations we 
not only use the impact parameter of the form but also D^, therefore 

the parameters for the two polarizations are given as: 

1 .|(l±^M.„)^B^^^^.O(^) ,4.94) 

4.2.3 General Trajectory to the Critical Orbit. 

Now that we have the general orbit equation (I4.92jl we can first test whether, 
for the critical impact perimeter D, the equation ^ — ^ for first order 
in A. Substituting the critical orbits and the impact parameters given in 
()4.54|) . ()4.55|) and ()4.94j) into ()4.92j) and expanding up to first order in A we 
have: For planar polarization. 

,du^o uf. , , , r / / AM ,,,, , AM ,,r, 

{j^f = ^(l + AMno)-[l-2MK(l + -^no))[(no(l + -^no)]2 

, ut,{l+\Muo)^^ ..^ ,^ AM AM 5 



-)[1 - 2M[no(l + --no)]][no(l + --^o)]^ 
o o 



(4.95) 
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and similarly for vertical polarization. 

(^)' = ^(1-AMmo)-[1-2M[mo(1--^Mo)]]K(1--^Mo)]' 

_ _o) _ 2[(mo(1 - -^Mo)]]K(l - -^wo)]' 

^ (4.96) 

Expanding these and eliminating all terms of order A"^ and higher, we find 
that the right hand sides go to zero, as required. Thus for the appropriate 
impact parameters these equations behave as they should. 

The next step is to solve equation ()4.92|) for the two polarizations. This 
is most simply done using numerical methods in Mathematica. As a guide, 
we know our solution will be of the form u = uq + kui, where Uq will be 
the classical solution ()2.39j) . ui will be a small modification that pushes the 
critical orbit up or down depending on photon polarization, and k will be 
some constant that is first order in A, i.e. of the form k = As, where s 
will be some number given by the boundary condition: for (j) ^ oo then 
u{(f)) no(l ± ^uo). If we substitute for the critical impact parameter D 
from ()4.54l) and ()4.55|) depending on the polarization, and then transform to 
u uq + Asui, and expand to first order in A, we have non-linear first order 
differential equations in uq, ui and 0. 

• For planar polarization: 

dui{(f)) . 1 



27M^uo{^y 



6^2717^ - ^o^*^)' + 2Mno(0)3 27M' 
+ 5AM\o{(pf - 6suo{(p)ui{(l)) + 18sMuo{(pfui{(j))] (4.97) 
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• For vertical polarization: 



d(j) 




54M^Uo{(pf - 6sno(0)ui(0) + lSsMuo{^fui{^)] (4.98) 



These can be solved in one of two ways, (i) is to substitute the classic so- 
lution ()2.39l) for uq and solve analytically, (ii) is to solve ^ and the classic 
equation for ^ simultaneously using numerical techniques. We attempted 
to use method (i) to derive an analytic solution for however, due to 

the complex nature of the equation we proceeded to use method (ii), that is 
solving by numerical methods. To do this we set the constant s = 1, and 
then when the numerical values of Mi(0) were determined we could deter- 
mine the constant s so that u{(f>) coincided with the modified circular orbits 
given in (|4.54|1 and ()4.55p . This technique was used for reasons of conve- 
nience, because solving (|4.97|] and (I4.98jl for various values of s would be 
time consuming as each numerical calculation takes a significant amount of 
time; therefore solving them once and then scaling the solution is a more 
convenient method. The results of the equations were plotted as: 



where the constant s was picked to satisfy the condition: — > oo =^ m(0) 



polarization. In this way the constant s was determined to be s = 1/3, 
which was tested for various values of M. We have plotted the results of the 
numerical calculation in Figs. 14.21 and 14.31 In Fig. 14. 2[ you can see that the 
general critical orbits follow a classic style path, however the orbit splitting 



m(0) = Mo(0) + AsUi{(f)) 



(4.99) 




Mo) i.e. the trajectories tend to the critical orbit, depending on 
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Figure 4.2: The quantum modified critical orbits follow classic type paths. 
We have used a factor A = 0.00009 in order to make the quantum corrections 
more visible, where a real value should be of order ~ 10~^^, Eqn. 14.381 The 
splitting of the orbits can be clearly seen in the Fig. 14. 3| given below. 



is not clearly visible. But in Fig. 14.31 we have plotted a closer view of the 
critical orbits, and here the splitting is highly visible. It can be seen that the 
general trajectories for the planar and vertically polarized photons tend to 
the relevant critical orbits. 
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^^^^^^^^^^^^^^^ FtenfifK]laraed(c<l> 










^^^^^^ '~ verlkal/ poldrt^ 



Figure 4.3: The quantum modified trajectories show clear sphtting as they 
tend to the critical points. We used the constants M = | and A = 0.00009. 

4.3 Quantum Modification and the Event Hori- 
zon 

In Chapter. H it was shown that when we decreased the impact parameter 
from the critical value {D = 3^/3M to D — 1/10) the photon trajectory, 
given as spiraled into the singularity, Fig. 12.51 and when we represent 

the trajectory as a function of coordinate time, u{t), it tends to the event 
horizon, u = 1/2M. From the horizon theorem it was seen that quantum 
modifications have no effect on photon velocities directed normal to the event 
horizon. So this implies when a photon tends to the event horizon at an an- 
gle, e.g. with an impact parameter D = DcnUcai — 1/10, then the component 
of velocity normal to the horizon should be unchanged, while the component 
parallel to it is modified according to the quantum correction; this modifica- 



CHAPTER 4. QUANTUM MODIFIED TRAJECTORIES 55 

tion should then result in a shift of the photon trajectory, but the horizon 
should remain fixed. In order to test this we used Eqn. (|4.92p with an im- 
pact parameter D = D critical ~ 1/10 to show that the photon trajectories 
still fall into the singularity. We then used Eqn. ()4.93|1 . with the new im- 
pact parameter, to study the behavior of the trajectories around the event 
horizon. 

4.3.1 Trajectories to the singularity 

In order to construct quantum modified trajectories, which go past the crit- 
ical orbit and fall into the singularity, we require the impact parameters: 



1 1 1 , 
— = y= ,D = —T- + 0{A) 4.100 

where, as before, + is for vertical polarization and — is planar polarization, in 
1/ D"^. For these impact parameters we numerically solved Eqn. ()4.92|1 . and in 
Fig. 14. 41 we can see that all the trajectories follow a classic type path into the 
singularity^. However, near the singularity you can see the splitting of the 
orbits as they tend to m — oo. In Fig. 14.51 we have shown a magnified view 
of the point where the trajectories cross the event horizon. In this figure 
it can be seen that the planar polarized photon (c < 1) crosses the event 
horizon at a point before the classic trajectory and the vertically polarized 
photon (c > 1) crosses it at a point after the classic trajectory. This makes 

sense, as the planar polarized photons are pushed towards the black hole 

^The quantum modifications to the classic trajectories are very small, and even if we 
use the hugely exaggerated value of ^ = 0.00009, as was used in Figs. 14. ^il and 14.21 the 
modification is hardly visible. So, in order to magnify the quantum correction even more 
we used A = 0.0009. 
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Figure 4.4: The quantum modified trajectories follow similar paths, as the 
classic trajectory, to the singularity, with M = | and A = 0.0009. 

and vertically polarized trajectories are pushed out, the vertically polarized 
ones must spiral further around the black hole to reach the event horizon 
compared to the classic or the planar polarized trajectories. 

4.3.2 Fixed Event Horizon 

To demonstrate the fact that the event horizon remains fixed atu = 1/2M we 
applied the impact parameters given by (|4.100|1 to Eqn. ()4.93j) . Again, solving 
numerically we found that the trajectories tend to the event horizon in the 
classic way. Fig. 14. 6| however, they are again slightly shifted. In Fig. 14. 7L 
a close up of the point where the trajectories tend to the event horizon at 
u = 1/2M, you can clearly see that the quantum modified orbits tend to the 
horizon before the classic orbit. This can be understood by the fact that the 
planar polarized trajectory is pushed towards the black hole, hence it has less 
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Figure 4.5: The classic and quantum modified trajectories crossing the event 
horizon u = 2M, with M = | and A = 0.0009. 

of a distance to propagate before it reaches the horizon, and even though the 
vertically polarized trajectory is pushed outwards, the fact that it has a faster 
velocity than c = 1 it reaches the event horizon before the classic trajectory. 
The other more important thing we can note from this calculation is that the 
event horizon remains fixed at the classic point u = 1/2M, due to the fact 
that, as was stated earlier, the quantum correction has no effect on photon 
momentum normal to the horizon. This means even though the quantum 
correction implies velocities greater than the speed of light, no photons can 
escape from the event horizon, so the black hole remains black. 
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Figure 4.6: The quantum modified trajectories follow similar paths, as the 
classic trajectory, to the event horizon, with M = | and A = 0.0009. 




t *' 



Figure 4.7: The classic and quantum modified trajectories reach the event 
horizon at diff'ering times, with M = | and A = 0.0009. 



Chapter 5 

Quantum Modified Schwarzschild 
Metric 

Now, using our previous results, we can construct a new metric that en- 
compasses the quantum corrections due to vacuum polarization. This metric 
should then be a sum of the classic Schwarzschild metric and the polarization 
dependent quantum corrections we derived in Eqns. ()4.88j] and ()4.89j] : 

G^^uk^k" = + = 

= (^^, - ^i?^,,AaV)A;^r = (5.1) 

5.1 Construction of the Metric 

Using the results in Sec. 14.2.^ we can write the first components of the quan- 
tum modified metric Q^^ and Q^^, for vertical and planar polarization respec- 



59 



CHAPTER 5. QUANTUM MODIFIED SCHWARZSCHILD METRIC 60 



tively, as: 



9tt 
F - 

F - 
9tt 

F - 
F - 



mi r r-^ 



AMD"' 



9 



-Fu^ 



-RtbtdO'2'^2 



nit 



8ca,2M DF^ 



.MF 



AMD^ 



)( 



9 



{3D'F\'> - Fu' 



(5.2) 



(5.3) 



where, in the last lines of each component, we have used Eqn (|4.31|) to replace 
the constants with A, and we've made the transformation r — > 1/u. Now, 
doing this for all the other components we can construct the metrics for 
vertical and planar polarizations: 



Vertical Polarization: = u{0, 0, 1, 0) 



/ F- BFu^ 



v 
















-4 





Planar Polarization: = n(0, DF, 0, D^) 
( F + BFu^d-SD^Fu"^) 









Q 









J_ _ Bu^{l~D^Fu^) 
' F F 



Bu^D^ 








(5.4) 



Bu^D^ 

dtp 



-4 - BD^Fu^ 
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where 5 = These are now the relevant quantum modified Schwarzschild 

metrics for vertical and planar polarizations^. 



5.2 Dynamics with the Quantum Modified Met- 
ric 



We can now derive the same equations of motion as ()4.92|) . but now we can 
do it simply with the quantum modified metrics. Using the general wave 
vector in a plane: 

where we have transformed to r ^ and we substitute for i and as 
before. Applying this wave vector to ()5.4|) a.nd ()5.5l) we find: 



= g'^^kv 

for vertical polarization, and 



= gf^^n" 

for planar polarization. These are identical to the equations we derived in 

secsra 



^In both the metric, G^^, and vector, ai^, we have made the substitution r ^ 1/u and 
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Radial Geodesies 

We can now show that the metric for vertical polarization is consistent with 
the fact that radially projected photon trajectories are not modified. By 
using the vertical polarization metric, ((531), ^ general radial wave vector, 
(A;^ = (t,f,0,0), we find: 



= 

dt 1 , , 



which is identical to the classic radial geodesic equation, ()2.3flj) . However, 
the same is not true for the planar polarization metric. Due to our deriva- 
tion of the polarization vectors in Sec. 14.2. IL we constructed a very general 
vertical polarization vector and then normalized it; however, the one for pla- 
nar polarization was constructed for the case where ^ 7^ 0, as can be seen 
in Eqns. M.BOjl and (I4.76jl : therefore the planar polarization vector has 
dependence "mixed" into it through the substitution of 0: 



(5.10) 



We can note that it's a result of the substitution ^ = ^, in the polarization 
vectors, that we acquire an extra parameter D in our quantum correction 
(apart from the one introduced through normalization, which was incorpo- 
rated into B = AD'^/9). We can then say that if we need to use the metric 
for a radial trajectory with ^ = we can just set D = 0, as this would 
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lead to the removal of the dependence. In this way, the planar polarized 
quantum modified metric for radial trajectories is: 



( F + BFu^ 






\ 

F F 





(5.11) 




-4 
/ 

Then applying this to a general radial wavevector, we again find the classic 
radial geodesic equation: 



= ^(0=o);.Fr 

dt _ 

dr ~ F 



(5.12) 



Therefore, we have metrics for the Schwarzschild space time that incorpo- 
rate quantum corrections due to vacuum polarization; and these metrics are 
consistent with classic results. 



Chapter 6 
Summary 



For Schwarzschild Spacetime we showed, in the orthonormal frame, that, 
due to quantum corrections, the stable circular orbit at uq = 1/3M splits 
depending on the polarization of the photon, and the new modified circular 
orbits are given by the classical orbit plus a correction term to first order in a 
dimensionless constant A: u = Uo+Aui, where A = '^^2°^2' and has an order of 
10"^^. As stated by the polarization sum rule this splitting of the critical orbit 
is equal in magnitude but opposite in sign for the two polarizations. We found 
that the vertically polarized photon (c>l) is pushed out by a correction of 
ui = — ^Uq, while the planar polarized photon is pulled in by a correction of 
ui = ^u^. It was also found that this orbit shift also requires an appropriate 

2 

modification of the classic impact parameter 7^ = This modification, 
for the quantum corrected orbits, took the form 752 = 752 + ADi, and again 
as for the orbit shift, the change was equal in magnitude but opposite in 
sign for the two polarizations: Di = for planar polarization and Di = 
_Muq vertical polarization. Using this information, for the splitting of 
circular orbits, we then constructed the quantum corrected general equations 
of motion for the Schwarzschild spacetime. Using these equation of motion 

64 



CHAPTER 6. SUMMARY 65 

it was shown that a photon starting at m = 0, with the appropriate critical 
impact parameter, tends to the critical orbit associated with that impact 
parameter - and the trajectory follows a similar path to the classic case. 

We then went on to show that, using the general quantum corrected equa- 
tions of motion, a photon projected towards the black hole with an impact 
parameter less than the critical value falls into the the singularity, in terms 
of the angular distance (0) . Although the photons follow a classic type path 
into the singularity, the trajectories are slightly shifted according to polar- 
ization. The planar polarized photon crosses the horizon before the classic 
trajectory, and the vertically polarized one crosses it after, this corresponds 
to the fact that planar polarized photons are pulled towards the black hole 
and vertically polarized ones are pushed away, hence the vertically polarized 
ones need to go a further angular distance to reach the event horizon. 

In terms of coordinate time (t) we found that the photon trajectories 
tend to the event horizon. Therefore, the quantum corrections do not shift 
the classic event horizon from u = which corresponds to the horizon 
theorem. Also, we found that, although the quantum corrected orbits follow 
a classic type path to the event horizon, the point at which they hit the 
horizon is, again, slightly shifted depending on polarization. However, this 
time both polarizations hit the horizon before the classic trajectory. The 
planar polarized photon tends to arrive at the horizon first, then the vertically 
polarized one, and finally the classic photon. This could correspond to the 
fact that the planar polarized photon, although it has a velocity lower than 
the speed of light, has less of a distance to go, as its trajectory is pulled 
towards the black hole. For the vertically polarized case, even though it has 
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a faster than light velocity, it has a further distance to go to reach the horizon 
as its trajectory is pushed away from the black hole. 

Having determined the equations of motion, with the quantum correction, 
we then went on to construct a Schwarzschild metric that incorporates the 
quantum correction: Q^^, = {g^^, + 0^,^), where the correction: g^^, was again 
first order in A. We showed that with this metric and a general photon 
wave vector, /c^, we obtain the quantum modified equations of motion, as 
before. Also, this new metric confirms the horizon theorem, that is, when 
we use a wave vector indicating a radially projected photon we obtain the 
classic equation of motions. This was fine for vertical polarization, however, 
in the planar polarization case we had a problem; we had previously used 
a substitution that mixed a "hidden" radial angle, </>, into our polarization 
vector (as in general orbits the planar polarization depends on 0). By tracing 
back to the origin of this substitution we found that, in order to study radial 
trajectories, we need to set D = 0, this then removes the (p dependency; 
this then gives us the classic equation of motion for planar polarized radially 
projected photons. 

So in conclusion, after studying the dynamics of null trajectories in Schwarzschild 
spacetime we derived the polarization dependent photon trajectories to first 
order in the constant A (which is dependent on the fine structure constant, 
the mass of the star, mass of the electron, and the energy of the photon). We 
then incorporated these modification into a general quantum modified metric, 
which could also be used to derive the general quantum modified equations 
of motion. Also, the results of this work coincide with the conditions of the 
horizon theorem and the polarization sum rule. 
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Superfluid Behaviour of the 2+ Id 
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Chapter 7 
Introduction 

7.1 Quantum Chromodynamics 

Since the 70's it has been accepted that nucleons and other hadrons (baryons 
and mesons) observed in particle accelerators are not fundamental parti- 
cles themselves, but are composed of fractionally charged fermions known as 
quarks, which exchange bosons of the strong force known as gluons. In this 
description, known as the quark model, baryons and mesons are depicted as 
bound states of three quarks and quark anti-quark pairs respectively. In this 
way the quark model provides a very natural explanation for the multiplicity 
and pattern of all the strongly-interacting particles [121 CHI. Experimental 
tests of this theory (in a similar way to the classic high-angle Rutherford scat- 
tering of a-particles off atoms demonstrating the existence of the nucleus) 
consists of high energy inelastic scattering experiments of electrons off nu- 
cleons. The evidence from such experiments is consistent with the presence 
of pointlike spin-| constituents called partons with a mass one third that of 
the proton. These partons, which are able to move freely within the nucleon 
volume, are then identified with quarks. 
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7.1.1 QCD: A Model of Strongly Interacting Matter 

Quantum Chromodynamics (QCD) was introduced in the early 70s as the 
theoretical framework that translated the experimental and conceptual de- 
scription of the quark model into a quantitative calculational scheme. QCD 
describes quarks and anti-quarks as quanta of the elementary fermion fields 
ijj and ip, each with an SU{3) colour charge, and gluons as quanta of a self- 
interacting non-abelian gauge field A^. The Lagrangian density of QCD is 
given by 

cqcd = m^jp - mo)rfi^^ - ^^;.-^r (7.1) 

In the fermionic part i and j run over Nf fiavours of quarks, a and (3 run 
over the 3 colours, then mo is an Nj x Nf mass matrix in fiavour space. The 
covariant derivative 

jp = = YiS'^^d, - '-g{XT^A^^) (7.2) 

is introduced so the Lagrangian density is invariant under local SU (3) gauge 
transformations, where g is the bare coupling constant, is a vector (gauge) 
field with eight gluonic degrees of freedom and A" denote the generators of 
the SU (3) group. Due to the introduction of the eight gluon fields, through 
gauge symmetry, we adjoin the free gluon lagrangian (the final (gauge) part 
of the action) to give the full QCD lagrangian. In the free gluon term of the 
action the field strength tensor is 

= d,A: - d^Al + gia,cA\Al (7.3) 
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where fate are the structure constants of the SU{3) group^. As the QCD la- 
grangian is symmetric under a non-Abelian gauge group, underlined by the 
presence of the structure constants /"*'^, the theory has some non-trivial fea- 
tures that are not present in Abelian gauge theories like quantum electrodynamics 

• '^QCD contains gluonic self-coupling (three and four gluon vertices), 
which means the gluons themselves carry colour. 

• At large momentum, Q, the QCD coupling behaves as: 

where Aqcd ~ 200MeV is the QCD scale parameter. 

• Eqn. (I7.4jl implies that as(Q^) ^ as — > oo, and asiQ"^) oo 
when = a^^^d- 

This behaviour of the strong force is known as asymptotic freedom, and can 
be simply represented through the quark anti-quark potential 

V{r) = + Kr (7.5) 

where K is an experimentally determined constant, called the string tension, 
with an estimated value of ^ f420Mel^)^[T3]: and is the coupling given in 
Eqn. ()7.4|1 . which varies with distance as 1/ ln(r~^). Now it can be seen that 
for small distances the first term dominates, in which case the strong force 
behaves like an attractive Coulomb potential, and in the limit r — the 
quarks can be considered as free non-interacting particles^^. However, with 



i[T', T"] = i/'™"r", where T'' are the generators of the SU{N) Lie group, 
^in units of energy over length we have K ~ 880MeV/ fm. 
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greater separations the potential scales approximately linearly due to the self 
interaction of the gluons, as seen through the second term. This is related 
to the phenomenon of "confinement", i.e. to the empirical fact that coloured 
objects, like quarks and gluons, do not exist as physical degrees of freedom in 
the vacuum. So, as the coupling constant (related to the potential) becomes 
larger for greater separations perturbative treatments of QCD become less 
and less effective. 

7.1.2 Lattice QCD 

A perturbative treatment of QCD leads to a successful description of the 
force between quarks at small distances . However, as already stated, at 
large distances a perturbative treatment of QCD is less fruitful. One method 
to address the non-perturbative nature of QCD at large distances is that of 
Lattice Gauge Theory, proposed by Wilson in 1974. In this method all the 
fields are defined on a discrete Euclidean space-time lattice with a nonzero 
lattice spacing a. Thus, in this way the lattice QCD calculations can be 
numerically carried out without the use of perturbative expansion. Even 
with its advantages, lattice field theory still has problem, such as: 

• It's difficult to discretise the fermion field in a chirally symmetric way. 

• Simulations with reasonable light current quark masses are computa- 
tionally very expensive. 

• To simulate a smooth space time by implementing a large lattice volume 
and small lattice spacing requires a very large number of lattice sites. 
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• Simulations for non-zero chemical potential (/i 7^ 0) are next to im- 
possible due to the sampling weight (used in Monte-Carlo methods) 
becoming complex fAppendix IC.ljl . 

However, with the continual improvement of lattice algorithms and the ad- 
vances in computing power, lattice QCD is the driving force in our under- 
standing of strongly interacting matter. 

7.2 Chiral Symmetry in QCD 

An important feature of QCD is its chiral symmetry for fermions with a 
vanishing mass (or approximate symmetry as is the case for physical quarks) . 
Chiral symmetry is related to the symmetries associated with a particle's 
handedness, which in turn is defined by its helicity. Before we go into chiral 
symmetry we'll briefly discuss the concept of helicity and chirality. 

7.2.1 Helicity and Chirality 

A particle propagating with spin s has helicity h = s ■ k/\k\, which is the 
projection of the spin axis along the direction of its motion k, where positive 
helicity is right-handed and negative is left-handed^. Also, if the particle has 
a vanishing mass its helicity would then be invariant under Lorentz transfor- 
mations, which (for the massless case) leads to two good quantum numbers 
Bl and Br, referring to the separate conservation of left- and right-handed 
particle numbers in the absence of external effects. However, in the case 

of massive particles these quantum numbers are not separately conserved, 

^So for spin ^ particles there are two possible helicity eigenstates given as: h — 
known as left- and right-handed states. 
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but their sum is a good quantum number: B = + B^. So, through this 
concept of helicity we can define the chirality operators that project out left- 
and right-handed field states^ 

= ^(1 - 75)^ = Pl-ip '^R = ^(1 + 75)^ = Pr-^P 

= V^(l + 75) =^i"R ^/;^ = ^/;^(l-75) =V5Pi (7.6) 
and these chiral field states, ipL and ipn, satisfy the equations 

75^L = -^L Ibi^R = +'<pR (7.7) 

where 75 = i7°7^7^7'^ and its eigenvalues are called "chirality"^. In general, 
unlike helicity, chirality is not directly measurable. We find that in the limit 
where m — > (or for E ^ m) the helicity and chirality of a particle are in 
one-to-one correspondence, and the chiral fields associated with the massless 
particles thus represent physical states. However, as helicity is not Lorentz 
invariant massive particle fields must be expressed as a sum of left- and 
right-handed chiral fields, 

4; = Pl^ + Pr^Ij = ^Ijl + ^r (7.8) 

which is the covariant formulation of massive particles. So in this case, of 
massive particles, chirality and helicity are distinct things, and thus the chiral 
states are not physical states, but represent internal degrees of freedom. 

7.2.2 Chiral Symmetry 

With this discussion of handedness and chirality we can go on to discuss what 

is meant by chiral symmetry. In general, chiral symmetry is the symmetry 

''The field operators which create and destroy a quark are ijj and respectively. 
^The gamma matrices also satisfy 7575 = 1 and {7^,7^} = 0. 
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associated with the independent transformations of the left- and right-handed 
chiral states of a particle. So, when we say QCD, for Nf quark flavours, 
possesses chiral symmetry under SU^Nf)^ x SU{Nf)ji^ we actually refer to 
the lagrangian of the theory being invariant for the separate transformations 
of the left and right-handed chiral fields [T9] 

V; -^VtP, V = VlPl + VrPr, 

^ ^^V, V = vIPr + V^Pl (7.9) 

where Vl, Vr G SU{Nf). This is equivalent to the symmetry under the 
SU{Nf)v X SU{Nf)A, with the transformations 

SU{Nf)v : ^ ^ e-^-^-V^ ip 7/;e^»^" (7.10) 

SU{Nf)A : ^ ^ e^^l^Taiia^ ^ _^ ^g-^75ra«?a 

where are the generators of fiavor SU{Nfy. Chiral symmetry would be 

exact in the limit of Nf massless flavours, but for non- vanishing mass (mg 7^ 

0) it is explicitly broken from SU{Nf)v x SU{Nf)A to S'f/(iV/)y^ as the 

mass term in Cqcd mixes the chiral states and However, though 

quarks have non-vanishing masses, chiral symmetry is still a useful concept 

for the up/down quarks {Nf = 2) as the masses are very small and could be 

considered negligible compared to the QCD scale parameter Aqcd (and to a 

lesser extent for the inclusion of the strange quark i.e. Nf = 3). Thus, as long 

as the masses are small compared to the relevant scale of the theory one may 

treat SU{Nf)A as an approximate symmetry, so that predictions based upon 

the assumptions of the symmetry should be reasonably close to the actual 
^For Nf = 2 the generators are the Pauh matrices. 

'''It can be shown that Cqcd is invariant under SU{Nf)Y, but not under SU{Nf)A- 
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results[T6]. As this is an approximate symmetry we could ask: what are the 
observable phenomena related to chiral symmetry being approximate rather 
than an exact symmetry? One of the most obvious effects is seen through the 
non-zero (but small) mass of the pions, i.e. the Goldstone bosons associated 
with the spontaneous breaking of chiral symmetry for Nf = 2 (which will be 
discussed in the next section). 

Total Chiral Symmetry Group of QCD 

It can be shown that the QCD Lagrangian is also invariant under phase 
transformations of the left- and right-handed quarks {ul and (11), which is 
the U{1)l X U{1)r symmetry. So the total chiral symmetry group of the 
QCD Lagrangian for Nf = 2 flavors can be written as 

U{2)lxU{2)r = SU{2)LxSU{2)RxU{l)LxU{l)n 
= SU{2)v X SU{2)a X U{l)v x U{1)a 

(7.11) 

where the axial and vector U{1) rotations can be written as 

U{l)v ■■ ip e-^V i^e'^ (7.12) 

U{1)a: ^-^ e'^"-^^ ip ipe'^'^ 

Having the full chiral symmetry group of QCD we can discuss the physical 
manifestations of theses symmetries in nature. The U{1)a symmetry is (said 
to be) explicitly broken due to quantum fluctuations [21], but the U{l)v sym- 
metry is responsible for baryon number conservation, and hence is labeled 
as U{1)b- The pure unitary transformation SU{2)v corresponds to isospin 
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conservation, where the axial transformation, SU{2)a, alters the parity that 
is associated with a state. This means the manifestation of SU{2)a in na- 
ture would require that each isospin multiplet be accompanied by a mirror 
multiplet that has opposite parity. However, as no such multiplets are ob- 
served in nature it is an accepted view that SU{2)a is a broken symmetry 
of QCD, where the resulting massless Goldstone bosons are associated with 
the pions. We can also write the conserved currents associated with each of 
the symmetries (as given by Klevansky [221): 



SU{2)v : 




= ip'^^T^il) =^ Isospin 


Uil)v : 




= ip'Jfj.ip => Baryonic 


SUi2)A : 


jk 


= ijj'y^'y5T''tp =^ Chiral 


U{1)a : 




= ip'Jij.'j^ip =^ Axial 



(7.13) 

We have used the notation I for the currents, this is due to the fact that in 
this thesis J is reserved for source terms, as is the convention in condensed 
matter physics. 

7.2.3 The QCD Vacuum and Spontaneous Symmetry 
Breaking 

As discussed, due to the smallness of the current quark masses, QCD is said to 
possess approximate chiral symmetry. However, in the world around us this 
symmetry is spontaneously broken due to dynamical mass generation, occur- 
ring when quarks interact with vacuum quark condensates. In this section we 
will discuss this phenomenon of spontaneous breaking of chiral symmetry as 
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it will help in building a conceptual picture of the QCD phase diagram, which 
will then lead to a greater understanding of the superconducting/superfluid 
phases of strongly-interacting particles at high baryon density. 

Dynamic Mass Generation 

In metallic superconductivity it is seen that a small electron-electron attrac- 
tion (due to effects of the surrounding lattice) leads to bosonic particles of 
bound electron pairs, known as Cooper pairs. These bosons then form a 
condensate in the ground state of the metal leading to superconductivity. In 
a similar way we can assert that the ground state of QCD, or vacuum, is 
unstable with respect to the formation of a quark condensate. 

Quark-anti-quark pairs are created in the vacuum^, near the surface of the 
Dirac the binding energy of the ipip pair exceeds the energy needed to 

excite the quark anti-quark pair. Once excited the strong attractive interac- 
tion between them causes the bound fermion pairs to condense, which leads 
to the creation of an energy gap. The resulting vacuum quark condensate is 
quantitatively characterized by a nonzero vacuum expectation value 

<ijij>=<o\tlJL^pR + ^RiJL\o>^o (7.14) 

This nonzero expectation value signals that the vacuum mixes the quark 
chiral states leading to them acquiring an effective mass. This can be seen 
conceptually if you consider that a left-handed quark^ propagating through a 

vacuum can be annihilated hy ipi, leaving ipn to create a right handed quark 

^These fermion pairs have zero total momentum and angular momentum, thus they 
contain net chiral charge i.e pairing left-handed quarks with the antiparticles of right- 
handed quarks 

^for E :s> m we can assume that ?7i ^ so its helicity and chirality will be equivalent. 
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chiral state with the same momentum. As this continually happens, with the 
quark traveling through the QCD vacuum, its chiral state will flip at a rate 
proportional to < ipip >, which implies it would propagate just as if it had a 
mass^° [T3] . 

Spontaneously Broken Chiral Symmetry 

Through this dynamically-generated mass, called the constituent mass S as 
opposed to the current mass m, chiral symmetry is spontaneously broken. 
This spontaneous breaking of symmetry, occurring due to QCD's own dy- 
namics, leads to massless Goldstone bosons. For the case of Nf = 2^^ we 
end up with three massless Goldstone Bosons, which are identified with the 
isospin triplet of relatively light mesons, the pions, tt^ and vr*^. These pions 
are light, but not massless, which (as previously discussed) is a consequence 
of the fact that chiral symmetry was initially not an exact symmetry as 
quarks have a nonzero, but small, mass to begin with. 

This generation of (constituent) rest mass, and the associated chirally 
broken phase, is portrayed in Fig. 17.11 In this energy-momentum diagram it 
can be seen that, due to the pairing of the quark anti-quark pairs, an energy 
gap of 2S is created between the highest (quark) and lowest (anti-quark) 
states. This gap then represents a rest mass that is far greater than the cur- 
rent mass i.e. S ^ m, which can be interpreted as a physical representation 

^"A naive, but some what helpful, picture is that of a spoon being dragged through 
honey, in such a case it would seem to have a greater apparent mass (or inertia) due to 
the viscous drag of the honevfT^. 

^^From Goldstones Theorem it can be shown that the breaking of the symmetry 
SU{2)v X SU{2)a (3+3 generators) to SU{2)y (3 generators) gives rise to 3 massless 
Goldstone bosons. 
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Figure 7.1: ^p^p pairing instability leading to a chirally broken phase of the 
vacuum. 

of a vacuum with broken chiral symmetry 

In Sec. 17.31 we will go onto discuss how a variation in thermodynamic 
conditions could lead to the restoration of chiral symmetry, and how ex- 
treme variation then leads to exotic strongly interacting matter i.e. colour 
superconductivity / superfluidity. 

7.3 The QCD Phase Diagram 

When it had become clear that hadrons were indeed a state of conflned 
quarks and gluons, it was suggested that they should become deconflned 

at high temperatures or densities when the hadrons strongly overlap and 

^^The vacuum with restored chiral symmetry also has a gap, equal to 2m, however as 
it is very small we talk about approximate chiral symmetry. 

^^This section is a condensed summary of 'The Phase Diagram of QCD' by Simon Hands 
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T 

Figure 7.2: The phase diagram of H2O (not to scale). 

lose their individuality [H]. In this picture, we then have two phases, the 
"hadronic phase" where quarks and gluons are confined, and the quark-gluon 
plasma where they are deconfined. Such phases of strongly interacting matter 
and the associated transitions between them, with external thermodynamic 
control parameters like temperature and density, make up a map of The QCD 
Phase Diagram. 

7.3.1 A Classic Phase Diagram of H2O 

A classical example of a phase diagram is that of H2O, Fig. 17.21 In this 
diagram the thermodynamic control parameters are temperature T and pres- 
sure P, and the different manifestations (or phases) of H2O are separated 



CHAPTER 7. INTRODUCTION 81 

into three regions of ice, water and steam, with the separating boundaries 
between them marked by the eqiiihbrium coexistence curves P{T). Thus, a 
first order phase transition such as melting or boihng is observed when mov- 
ing along a path in the (T, P) plane that crosses the equilibrium boundaries. 
In this diagram there are two special points: the critical point (Tc = 650i^, 
Pc = 2.21 X lO'^Nm^'^) where the phases of water and vapour become indis- 
tinguishable, and the triple point {Ttr = 273.16K ,Ptr = 6007Vm~^) where all 
three phases coexist. 

7.3.2 A Map of the QCD Phase Diagram 

A possible schematic of the QCD phase diagram is presented in Fig. l7.3[ with 
the thermodynamic control parameters temperature, T, and baryon chemi- 
cal potential, /i. In this diagram we can note several locations which have 
analogues in condensed matter physics. The first location is the bottom left 
hand corner of the phase diagram where T and fi are both small. Here the 
thermodynamic behavior of QCD can be described as that of a vapour of 
hadrons, i.e. composite states of quarks and/or anti-quarks. In this loca- 
tion extensive work has been carried out to classify and quantify the bound 
states of strongly interacting matter, hence, this area has a certain amount 
of resemblance to relativistic atomic physics. It is then argued, and has been 
shown through numerical calculations^^ that as we increase T the hadronic 
vapour phase cannot persist. Eventually there comes a point, found to be 
Tc ~ 170Mel^[25]. where the dominant degrees of freedom are no longer 

hadrons, but quarks and gluons themselves. This "gas" of strongly interact- 
^''For a discussion of this see and references therein. 



CHAPTER 7. INTRODUCTION 



82 



Early universe 

Quark-gluon plasma 



Crossover 



n,= 
Vacoium 




Quark matter 



<\j; i|;> > 
Hadronic fluid 



Crossover 

Superfluid/Superconducting Phases 



" \ 2SC <HJV>=-0 CI L 

Nuclear matter i 

Neuttort star cores 



/ig- 922 MfV 



Figure 7.3: The proposed phase diagram of QCD (not to scale). 



ing matter is known as the quark-gluon plasma (QGP), thus the upper left 
region of the phase diagram has the classic analogue of relativistic plasma 
physics. Two possible candidates for the formation of the QGP phase are the 
early universe, immediately after the Big Bang, and high energy collisions 
between nuclear particles. 

The final area of the phase diagram, and also the most elusive, lies along 
the /i-axis. Unlike simulations along the T-axis, lattice gauge theory simula- 
tion become ineffective when applied to QCD with /i 7^ 0. However, through 
extrapolation methods it is seen that as you increase /i, for zero tempera- 
ture, the ground state (i.e. the vacuum) persists until /i reaches the value of 
nuclear matter (nucleon rest mass minus the binding energy), at which point 
it becomes energetically favorable to fill the ground state with a bound nu- 
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cleonic fluid. This onset value is estimated at /io = 922MeV, at which point 
baryon density Ub jumps from zero to nuclear density Ubo — 0.16/m~^. This 
is referred to in [THj as 'room chemical potential' as the vacuum and nuclear 
matter coexist at this point. Extrapolation methods, can be employed for 
densities up to 2 — 3^^0 5 however beyond this we are forced to rely on approx- 
imate treatments such as the MIT bag model or the NJL model. Through 
such treatments it is estimated that as increases we again expect a tran- 
sition from a phase where matter exists in the form of nucleons to one where 
the dominant degrees of freedom are quarks and gluons. This phase of dense, 
strongly interacting, matter is believed to exist at the cores of neutron stars. 
It is also speculated that under such conditions a phenomenon similar to 
the Bardeen-Cooper-Schrieffer (BCS) instability, that leads to superconduc- 
tivity in metals and superfluidity in liquid '^He at low temperature, results 
in the formation of diquark pairs and the onset of a colour superconduct- 
ing/superfluid phase (see Sec. 17. 4j) . Therefore, the lower right region of the 
diagram has similarities to a branch of condensed matter physics; and it is 
this region of the phase diagram with which this thesis will be dealing. 

7.3.3 Restoration of Chiral Symmetry 

Before we go on to coloured BCS phenomenon (the area with which this 
thesis is concerned) we will discuss the restoration of chiral symmetry in the 
various areas of the phase diagram, in particular the areas of high T and /i, 
i.e. the upper part around to the bottom right of the phase diagram. 



CHAPTER 7. INTRODUCTION 84 

< ipip >~ For Extreme T 

It has been seen through numerical simulations that the onset of QGP co- 
incides with the restoration of chiral svmmetrv|25|. As you move from the 
phase of bound strongly interacting matter, at low temperature and density, 
up to the high temperature region dominated by QGP there is seen to be a 
drop in the vacuum quark condensate i.e. the order parameter < ipip > tends 
to zero in the QGP phase. However, to be precise it has been shown 
that the order parameter < ipip > doesn't strictly vanish, but rather drops 
very steeply in the transition region, and this transition is seen to occur at 
Tc ^ nOMeV. Thus, in the QCD phase diagram the formation of QGP 
along the fi = axis is naturally referred to as a crossover rather than a 
true first or second order phase transition. In a hand waving way we can 
interpret the "destruction" of the vacuum quark condensate as a result of 
the high temperature infiuence on the strong coupling constant, g{T). We 
can say that at very high temperatures, where QGP occurs, large energies 
are exchanged in inter-particle collisions, which results in the weakening of 
the strong interaction due to asymptotic freedom. Therefore, as the energy 
needed to excite quark anti-quark pairs is no longer below the binding energy 
between them the quark condensate does not form, which in turn leads to 
the restoration of chiral symmetry. 

< ijj^ >~ For Extreme /i 

As previously discussed, to study the high density region of the QCD phase 
diagram, i.e. above 2 — Sra^o (and T ^ 0), one must employ alternative mod- 
els to lattice gauge theory such as the MIT bag model and the NJL model. 
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Each model has its pros and cons: the most important characteristic of the 
bag model is its property of quark confinement, and that of the NJL model 
is chiral symmetry and its spontaneous breakdown in the vacuum; on the 
other hand the bag model violates chiral symmetry and the NJL model does 
not possess confinement. So one must decide which property is of greater 
importance and employ the appropriate model. Through simulations, using 
the NJL model, it is seen that as you increase baryon density beyond the 
nuclear density n^o, i-e. move from the hadronic phase to that of strongly 
interacting matter, a phase transition takes place from the chirally broken 
phase to one of restored chiral symmetry ^2]. This restoration of chiral 
symmetry in the high density regime can be considered due to the quarks 
becoming the dominant degrees of freedom at this density. At extreme den- 
sities (/i ~ 1200Mel^ [ini) hadrons strongly overlap, lose their individuality 
and start to behave like one big mass of quark matter; which is thought to 
occur in neutron stars. Being fermions, the quarks occupy the Fermi sea up 
to the level Ep. This leads to ipip vacuum excitations, for \k\ <^ kp, becoming 
Pauli-blocked. Therefore, at some point, these pairs require so much energy 
to excite that it becomes preferable to revert to a chirally symmetric ground 
state. In this way, at high density, chiral symmetry is restored. 

7.4 Coloured BCS Phenomenon 
7.4.1 BCS Processes 

In BCS theory the pairing of Fermions (known as Cooper pairs) leads to 
the formation of Bosonic particles, which condense and leave an energy gap 
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Figure 7.4: ip^p pairing instability leading to super-conductivity. 



2A that is equal to the energy of the binding between the pair, Fig. 17.41 
In metals this BCS condensation of electrons leads to superconductivity, i.e. 
the flow of electric current without resistance. It is also seen that BCS con- 
densation is accompanied by the breaking of the U{1) local electromagnetic 
gauge symmetry. Such BCS processes are also relevant in the study of super- 
fluidity (the flow of fluid with negligible viscosity). One type of superfluid, 
'^He, is a Bosonic fluid, and simply undergoes Bose-Einstein condensation at 
low temperatures. However, ^He, a Fermionic fluid, undergoes a BCS type 
process which leads to the condensation of bound Cooper pairs. As both 
superconductor and superfluid have an energy gap, the way to distinguish 
between the two is that a superfluid is characterised by a ground state which 
does not respect a global symmetry of the underlying action. 
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7.4.2 BCS in QCD 

Due to the strong attraction between quarks the BCS process leads to some 
interesting results in QCD. One of these, as we have already seen in sec. I7.2.3[ 
is the condensation of quark-antiquark pairs. This leads to spontaneous Chi- 
ral symmetry breaking, and a resulting effective quark mass which is much 
higher than the current mass. Also, it is believed that at high densities, 
when chiral symmetry is restored, and quarks are deconfined, condensation 
of diquark pairs occurs; and the resulting qq wave function is a gauge non- 
singlet. This condensation leads to the formation of an energy gap separating 
the ground and excited states by 2A. As in "orthodox" BCS processes, this 
condensation of diquark pairs is accompanied by the breaking of local colour 
gauge symmetry, SU{3)c; and in analogy with electromagnetic superconduc- 
tivity this high density phenomenon is known as colour superconductivity 
(CSC). 

If we were able to conduct Lattice QCD simulations at 7^ it could 
be possible that we'd see the colour superconducting phase. However, due 
to the sampling weight, used in Monte-Carlo methods, becoming complex, 
to study quark matter in the high density phase we require effective field 
theories, such as the NJL model. The NJL model is a purely Fermionic 
model, containing no gauge degrees of freedom, and interactions are simply 
represented by a four-point interaction. One of the most important features 
that makes the NJL model an ideal effective theory of QCD at high den- 
sity is that it not only observes the same (Chiral) symmetries of QCD, but 
the breaking of these symmetries in the vacuum is analogous to BCS super- 
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conductivity. For this reason it is viewed as an appropriate model to study 
coloured BCS phenomena. However, in the NJL model the qq wave function 
is gauge singlet, which implies the ground state is not superconducting, but 
rather superfluid^^. 



^^Another method to simulate the high density region of QCD is using 2 colour QCD, 
which also has a diquark condensate gauge singlet state. 



Chapter 8 

The Nambu-Jona-Lasino Model 



Quantum chromodynamics is the accepted theory of strong interactions and 
at short distances is highly successful. However, as discussed in the previ- 
ous chapter, for larger distances where perturbative techniques break down 
lattice gauge theory is required to calculate the QCD expectation values. 
This method itself has problems such as the requirement of huge computing 
power, problems associated with including fermions on a discretized lattice, 
problems for fi 0, etc (See Sec. I7.1.2|) . For this reason it is reasonable 
to isolate the relevant physics associated with a process of interest and con- 
struct an approximate model of the exact theory that accentuates the main 
features of the theory; to study strongly interacting matter at high density 
(i.e. /i 7^ 0) we do just this. In this chapter we start off by discussing the 
Nambu-Jona-Lasino (NJL) model as a model of QCD, its chiral symmetry 
and other associated technical issues, and also the benefits of its use for lattice 
calculations at high density. We then go onto discuss the lattice transcription 
of the NJL model in 2 + 1 dimensions. Finally we discuss QCD calculations 
involving the lattice NJL model in 2+1 dimensions, referring to the results 
of [12] and how they relate to the investigations of this thesis. 

89 
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8.1 The NJL Model of QCD 

Non-perturbative QCD calculations in the high density regime are next to 
impossible due to the measure of the Euclidean path integral becoming com- 
plex for baryon chemical potential jj ^ 0. This phenomenon, known as the 
sign problem, involves the QCD path integral effectively generating mean- 
ingless complex probabilities. One way to overcome this problem is by the 
use of the pre-QCD model of Nambu and Jona-Lasino [23 Ell- this model 
nucleons were considered to interact through a two body interaction, and 
the resulting lagrangian was symmetric under the full chiral symmetry group 
given in Eqn. ((ZTTJ. However, in accordance with experimental observation 
the U{1)a symmetry was later excluded. 

As this theory was conceived before the advent of quarks, Nambu and 
Jona-Lasino used nucleons as the elementary building blocks of the model. 
In analogy to the effective electron-electron interaction in BCS theory the 
nucleon-antinucleon attractive interaction was considered to be responsible 
for the formation of Cooper pairs of nucleons and antinucleons, which results 
in a nucleon condensate. This condensate would then lead to a vacuum 
mass gap and the breaking of chiral symmetry, where the pion could then be 
identified as the Goldstone boson occurring in the theory due to the breaking 
of the axial symmetry. 

With the possibility of quarks the NJL model can be reinterpreted as 
a simple model of QCD, with the quarks taking the place of the nucleons. 
However in this simplified model the strong interaction, mediated by the 
gluons, is excluded; instead the interaction between quarks is assumed to 
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be a point-like four-point interaction (for Nf = 2). This means it does not 
describe quark confinement - the phenomenon whereby single isolated quarks 
are never observed. On the other hand the NJL model has other interesting 
properties which make it a suitable substitute for QCD, i.e. the model's 
global symmetries and the patterns of their breaking are similar to those of 

QCD [221. 

8.1.1 NJL Lagrangian 

In Euclidean space, the original lagrangian density describing the NJL model 
with Nf = 2 quark flavours in the isospin representation of SU{2) can be 
written as 



where ifj and ip are independent Grassmann Dirac 4-spinors representing the 
quark fields, mo is the current, or bare-mass of the quark, r = (ri, r2, Ts) is a 
vector of the 2x2 Pauli matrices, which run over internal isospin or flavour 
degrees of freedom, and Tq is the unit matrix. 

When the NJL model was reinterpreted as a quark model, with 'ip as the 
quark field with two flavour and 3 colour degrees of freedom, it was kept in 
the original NJL form (|8.1L and in this thesis we will work with this form of 
the model. However, Eqn. (18. 1|) is not unique and other chirally symmetric 
interaction terms can be included. A more general form of the NJL model is 
given as [TT] : 




(8.1) 



C = ^(^+mo)ro^>^-|-[(V^ro^V)'-(/^''^'^)'+(/75ro^V)' 



(8.2) 
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8.1.2 Bosonization 

Through the introduction of auxiliary scalar and pseudo-scalar fields, a and 
TT, the NJL model can be rewritten in a bosonized form, which makes it easier 
to treat both numerically and analytically. Starting with the NJL Lagrangian 
()8.H) we can introduce the generating functional 

Z j VipViPexp[j d^xCiip,^;)] (8.3) 

Then using the fact that path integrals of Gaussian functions can be per- 
formed exactly j22|, 

J V^exp[ J d'^x{±A^ - B^^)] ~ exp[ J d^x^] (8.4) 

we can write: 

exp[ j d^xG{i)^f] ~ j PcTexp{ j d'^x[{i)i))a - ^f^^]} 

explj d'^xG{tjji'j5f%jjY] ~ J VirexplJ d'^xKtpi'j^f ■ tt^Jj) — ^^n'^]} 

(8.5) 

where we have G = — y- Therefore, we can rewrite ()8.3p in a bosonized form 
as 

Z ~ y VtjjVijjVaVTf exp{ J d'^x[%l)[(^ + mo + (o" + ■ 

- 4^(^' + v?-vf)]} (8.6) 

Now using $ = 0" + ivf ■ r and cr^ + vr ■ vr = |Tr$1^$ we have the bosonized 
form of the NJL lagrangian: 

C = i}[^ + mQ + {a + i-f5T ■ vf)]^ - -^rr$"^$ (8.7) 

oG 
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8.1.3 Chiral Symmetry of the NJL Model 

The simplest way to show that the NJL lagrangian is symmetric under chi- 
ral transformations ()7.9|) is to use the bosonized lagrangian in ()8.7j] . For 
vanishing mass the kinetic term can be easily shown to be symmetric^: 

^^d^^ ^ ijVYd^V^ = ^r{VlVLPL + VlVRPR)d^iP 

= i'Yd^^ (8.8) 

However, the symmetry of the bosonic terms is a little more subtle. Bosonic 
fields, as in the last term of (|8.7|) . transform as|19|: 

$ ^ Vl^vI (8.9) 

Using this we can show that the first bosonic term is symmetric, by first 
writing it in the form: 

V'(a + Z75f ■ 7f)V^ = ^l)L{cf + IT ■ Ti)%l)R + ^ii{a - if ■ T;)ipL 

= ^L^ijR + ^R^^ijL (8.10) 

where 

\ —712 + ^TTi CT — ZTTs J ^ ' 

Now under the chiral transformation we have: 

^jyl^VR^R + V'i?V^$ViZ^i (8.12) 

which is symmetric, as the bosonic field $ transforms as (I8.9jl . 

^VV ~ V^VlPl + vIVuPr implies mass term breaks chiral symmetry fSec. 17.2.211 . 
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8.2 The Lattice Transcription of the N JL Model 



In this thesis we use the lattice transcription of the NJL model in 2 + 1 
dimensions, as used in [l2]. We replace the space-time continuum with a 
3-dimensional lattice, where each site is separated from its nearest neigh- 
bours by an arbitrary lattice spacing a. In dimensions d > 2 + 1 this lattice 
spacing must be chosen for the effect of introducing an ultra-violet (UV) 
cutoff A ~ a~^, which then regularises the divergences and makes the theory 
mathematically well-defined, however, in ci = 2 + 1 dimensions the theory is 
renormalizable and well defined on the lattice for an arbitrary cutofffT^^. 

The structure of our model consists of the fermion fields being defined on 
the lattice sites x and the bosonic fields being defined on the dual lattice sites 
X, translated from the original lattice by (|, |, |). The space-time integrals 
and differentials are represented by sums and finite differences, so that in the 
staggered fermion notation the noninteracting quark action 



where x ^ind x ^ire the Grassmann-valued staggered fermion fields defined 
on the lattice sites, and r]^{x) is the Kawamoto-Smit phase i)^o+ 
Now, using this definition of the free action, we can go onto describe the 

construction of the full NJL model as used in this thesis. 

^As the interaction between quarks is assumed to be a point-like in character, the 
theory is not renormalizable for dimensions greater than 3 + 1 . 




(8.13) 



is given by 




(8.14) 
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Lattice NJL Model for /i = 

The NJL model on the lattice for /i = 0, in its bosonised form, is defined by 
the Euclidean action 

X ^ X 

where x, X ^^re, as before, the fermionic fields defined on the lattice sites, and 
$ = cr + il^ ■ is a 2 X 2 matrix of bosonic auxiliary fields, again as before, 
defined on the dual lattice sites. Also, M[$] is the kinetic operator, and is 
in the form for staggered lattice fermions interacting with scalar fields, 

1 ^ 

u=0 

<x,x> 

Here the parameters are bare fermion mass m and coupling constant g^. The 
Pauli matrices, r, are normalized to tr(rjrj) = 25^ and act on the internal 
SU{2) isospin indices p,q = 1,2, and e{x) denotes the phase (— l)^o+a;i+x2_ 
Finally we have < x, x >, which represents the sum over the set of 8 dual 
lattice sites neighbouring x. In this bosonized form we can integrate over the 
auxiliary $ fields which leads to an equivalent action in terms of fermions 
that self-interact via a four-point contact term proportional to g"^, as shown 
in Sec. [HX21 

Lattice NJL Model for ^ 

To introduce a chemical potential /i into a fermionic field theory we must 
incorporate a term of the form £^ = fiip'jQip into the fermionic Lagrangian. 
The inclusion of this term can be understood in a simple way by considering 
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that to introduce a source into a Lagrangian we must incorporate a source 
term of the type Csource = J(p, where J is the source for the field ip. Therefore, 
as ip'joip ~ ip'^ip is the quark number operator Ng, we can think of as a 
source for this operator, which leads to the creation of fermions over anti- 
fermions as a more energetically favourable process. Therefore, the chemical 
potential term would then be incorporated into the partition function in 
the usual way, as e"^^'. However, introducing a chemical potential in a 
discretized field theory is not as simple. We cannot simply implement a term 
of the form given above and naively discretise it, as it can be shown that this 
leads to divergences in the energy density of the theory [23 EH]. This can be 
resolved if we consider that incorporating a chemical potential in the form 
given above is like the zeroth component of a vector potential interaction 
term in a fermionic Lagrangian, tp'jf^ipAf^ =^ ip'jQipA^. Thus, due to this 
similarity to a gauge field, we must incorporate the chemical potential term 
in a "gauge-invariant" way. This is done by introducing /i into the zeroth 
derivative by multiplying the forward time difference by and the backward 
time difference by e~^. In this way the kinetic operator ()8.15|) becomes 

+ l^xy Yl [^(^)^"' + '^(^M^) ■ (8-16) 

<x,x> 

Introduction of the Diquark Sources 

Due to limitations of computer memory we are forced to work on a finite 
volume system, and working under such conditions leads to certain compli- 
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cations. One such issue is that due to finite volume eff'ects on spontaneous 
symmetry breaking. A simple example of this phenomenon is seen in the 
0(N) model. One of the most important properties of this model is the 
effect of spontaneous spin symmetry breaking at the critical temperature, 
Appendix IbI In the 2-dimensional 0(N) model, above the critical temper- 
ature, the spins at each lattice site are randomly orientated due to thermal 
fluctuations, hence we have the average value < Si >= 0. When we hit the 
critical temperature, thermal fluctuations become less dominant compared 
to spin interactions and the correlation length increases, this in turn leads 
to spontaneous magnetisation (i.e. magnetisation with a zero external held). 
Then, at T = the magnetisation is at its maximum, and the correlation 
length is comparable to the size of the lattice. However, it is seen that such 
effects only occur in the thermodynamic limit 



These flnite volume effects are shown in Fig. 18. H which shows magnetisation 
verses external magnetic held. It can be seen that spontaneous symmetry 
breaking only occurs in the limit of ^ oo. So how can we simulate an 
inflnite lattice? This problem can be overcome in computational calculations 
by including a symmetry breaking term that represents the external field. 
Therefore, for a finite volume system we can conduct simulations with various 
values of the external field, such as the magnetic field H in the 0(N) model, 
and then extrapolate to the limit of the field going to zero, H ^ 0. This 
then gives the spontaneous magnetisation value of the system at temperature 
T. A similar problem is seen in the study of chiral symmetry breaking. 




(8.17) 
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where we include a finite bare mass m to allow the measurement of the 
chiral condensate < ipi/j > on a finite volume. It is for a similar reason 
we include the diquark and anti-diquark symmetry breaking terms, which 
allows us to measure the values for < ipip > and < iptp > condensates. 
The implementation of such diquark symmetry breaking terms requires the 
introduction of diquark and anti-diquark fields j and j that can be taken to 
zero for the infinite volume system. Therefore, in analogy with the mass term 
(for the measurement of the chiral condensate), and the external field term 
(in the 0(N) model for the measurement of spontaneous magnetisation) we 
include the following terms into the NJL Lagrangian: 

jyVaX + jM*'' (8.18) 
Final Form of the Model and its Symmetries 

We can now write down the full NJL model, as studied in this thesis: 

SnJL = Sfer + Sbos (8.19) 
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where the fermionic and bosonic parts of the action are: 

Sfer = J2 XM^X + JX*V2X + JXr2X''' (8.20) 

X 

Sbos = \y^tr<l>^^ (8.21) 

^ X 

where M[$] is given in Eqn. (|8.16jl . To construct the partition function of 
this system we rewrite the fermion part of the action by defining a bispinor 
= (x*^,x). In this, Gor'kov, basis we can write the fermion action as 

Sfer = ^''A^ (8.22) 

where the Nambu-Gor'kov matrix A is 

Apq — [ J '2 "xy 2^^^xy \ (R9^\ 

In this basis we can integrate out the fermion fields, which leaves us with the 
following Euclidean path integral: 

exp{—Sbosi^)) (8.24) 

The discretized NJL model, given by Eqns. dOn - lUTTT]) and (IOT)ll . still has 
the full QCD symmetry SU{2)l x SU{2)r xU{1)b, as discussed in sec. 18.1.31 
But now the projection operators (I7.6jl become Pr/l Pe/o = |[1 
which project onto even and odd sublattices, respectively; and the transfor- 
mations take the form: 

X ->(KPe + KPo)X X^X{PeV} + PoV^,) 

$ ^Vo^V} [Ve,VoeSU{2)] (8.25) 
X ^e^° X^Xe-^° [e*"et/(l)B] (8.26) 
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The symmetry of the diquark source terms, Eqn. ()8.18|) . requires an extra ex- 
planation. These terms were added to break baryon symmetry, U{1)b, hence 
allowing the measurement of the diquark condensate as j 0. Therefore, 
under the transformation ()8.26|) it is not symmetric, 

JX^V^x + JXr2r - JX*V2xe^'" + JXr2r e"^'", (8.27) 

where the symmetry is restored for vanishing j and j. Under the SU{2) sym- 
metry, ()8.25|) . the terms are still symmetric. We can show this by considering 
only the j term: 

jyV2X ^ JX*^(K*^Pe + CPo)r2(KPe + KPo)X 

(8.28) 

Now using the identity: T2UT2 = U* we can write^: 

Jx'''{r2VXPe + r2VXPo)x = 3X*'r2{Pe + Po)x 

= Jx''r2X (8.29) 

which shows the symmetry of the diquark source term (similarly for the 
anti-diquark source term). 

8.3 QCD Dynamics using the NJL Model 

Previous studies of the 2 + Id NJL model have shown that a strong first 

order phase transition, at yUc — S ~ 0.65, takes place to a chirally symmetric 
^Transpose of this identity gives: (— T2)C/*''(— T2) — U'' =^ C/*''t2 ~ T2UK 
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state [12]. This transition is seen to be accompanied by the onset of a non- 
vanishing density of baryon charge in the ground state. However, in this 
high density region it is also seen, for 2 + 1 dimensions, that there is no 
diquark condensation and the U{1)b symmetry is not broken, and there's 
no evidence of a non-zero energy gap. More interestingly, it is seen that the 
diquark condensate scales non-analytically at high density, < ipip >cx . 

In this section we will look at how these results were determined, and more 
importantly (for our work) how the non-analytic behaviour of the diquark 
condensate could imply superfluidity. 

8.3.1 Phase Transitions at High Density 
Chiral Condensate 

Chiral symmetry restoration is signalled by the first order transition of the 
chiral condensate, < ipip 0, which is determined as: 



1 d\nZ 
V dm 



1 dZ 
ZV dm 




(8.30) 



using det 2A 



exp(tr ln(2^)), we have: 




(8.31) 
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From (l8T6ll and ^^2^ we have 



dA f is 



dm V -^5 



.32) 



therefore (j8.31j) becomes 



(8.33) 



Then, the expectation we calculate is: 

< >=< ^tr^-i ( |^))> (8.34) 

Baryon Charge Density 

As chiral symmetry is restored at high density, and nucleons dissociate into 
quarks, there is seen to be an onset of a non-vanishing density of baryon 
charge in the ground state, ub =< V'ToV^ >, which is determined using Eqn. 

1 d\nZ 1 dZ 

riB = < tpioV > 



V dji ZV dji 
L j P$tp$ltr(^-i^)v/dit27exp(-5,,,($)) 

(8.35) 



where, as before, we find using ()8.16|) and (|8.23|) 



dA Q 



1 dhl 

2 dti 



(8.36) 
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where 

dM 1 

dM*'' 1 



Now, defining: 
X = ^ 



U 5P'[e''5y,x+o + e-%,x^d] 



4 V -S^ne%,x-o + e-^^,x+6] 



we then have: 

< i^'jo^ >-- 



.37) 



5.38) 



^/ ^^'^•^tr(^-iA') Vdet 2Aexp{-Stosm (8.39) 

We then determine the baryon density by the expectation: 

riB =< -^tiiA-'X) > (8.40) 

These transitions in the chiral condensate and the baryon density can be 
seen in Fig. 18. 2| from which we can note the strong first order transition at 
/x = 0.65. 

8.3.2 Diquark Condensation 

The diquark condensates < ipip > are calculated in a similar way to the 
previous expectation values, but this time the derivatives of the partition 
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Figure 8.2: Chiral condensate S and baryon charge density ub as a function 
of /i (from [12]). 



function are taken with respect to the diquark source (or the anti-diquark 
source for < >): 



< ilJip >-- 



1 d\nZ 



(8.41) 



V dj 

In order to effectively study diquark condensation the operators V'V'i are 
defined as: 



1 



< i^^± >= TTTT ( 



dhiZ _j_ dlnZ 



) 



42) 



2V ^ 

with corresponding source strengths j± = j ±j, so when j = j ^ j- = 0. 
Then diquark condensation is given by: 



< QQ+ >-- 



1 dlnZ 



< trro^"^ > 



^.43) 



V dj+ W 

Then, the non-vanishing of ()8.43l) in the hmit j+ ^ is a criterion for the 
formation of a diquark condensate and the resulting spontaneous breakdown 
of the f/(l)B symmetry. 
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Criticality and Superfluidity 

However, in p2j it was seen that the diquark condensate does not extrapolate 
to a non-zero vahie for high density, but scales non-analytically to zero as 

KifjifjXxj^, (8.44) 

implying no breakdown of U{1)b symmetry. It was then suggested that this 
behaviour could be due to working in 2H-1 dimensions. The Mermin-Wagner 
theorem predicts: there can be no phase transition to a state with long-range 
order at finite temperature for a two dimensional system with a continuous 
global symmetry [21]. This was supported in [111, where diquark condensation 
and an energy gap were observed in simulations of the 3+ld NJL model. 

It is also seen that the exponent 6 varies continuously with chemical po- 
tential, /i, taking the value 5 ~ 3 at /i = 0.8 and ~ 5 at ;U = 0.9. This 
suggests a line of critical points for /i > /ic. This type of behaviour is simi- 
lar to that observed in other 2-dimensional critical thermodynamic systems. 
For a spin system (XY model) at its critical temperature the spontaneous 
magnetisation, M, scales with applied magnetic field, if, as M oc , where 
again the exponent forms a line of critical points.^ 

The simplest description of a critical system is given by the XY model, 
for which long range order would also spontaneously break a U{1) global 
symmetry. Using this system Kosterlitz and Thouless (KT) described the 
critical behaviour, of a 2 dimensional system, as a consequence of binding 

and unbinding of vortex pairs at a critical temperature T^t (Appendix iBjl . 

''A similar non-analytic behaviour is seen for a fermionic model exhibiting chiral sym- 
metry breaking: < tpTp >oc m'. 
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Applying this to a 2-dimensional superfluid, the critical temperature is seen 
to be the point where superfluid vortices bind and unbind. Thus, above 
the critical temperature when the vortices are free they experience long- 
rang mutual interactions leading to the behaviour of a viscous fluid; but 
below the critical temperature vortex-antivortex pairs bind, reducing long 
range interaction, resulting in superfluid behaviour. The jump from the spin 
XY model to superfluids comes from the fact that we replace the total spin 
angle (of vortex-antivortex pairs) in the XY model by the phase angle of the 
possible superfluid condensate wave function (Appendix IB. 3|1 . Therefore, (as 
was conjectured in [12]) in analogy to the non-analytic behaviour displayed 
by the XY model the high density (and low temperature) critical phase of the 
2+ld NJL model can be considered to be a gapless thin fllm BCS superfluid. 

In the rest of this thesis we will investigate this possible superfluid be- 
haviour of 2+ld NJL model. This will be done by extending the baryon 
density to a baryon three current by the use of Ward Identities. Then, by 
the use of a spatially varying diquark source we will introduce a gradient 
in the diquark pair wave function, which should force a flow of the baryon 
current, and hence result in a measurable signal. Using this current signal 
we will explore its behaviour (using the helicity modulus, T, Sec. 19.21) with 
variations in spatial volume, temperature, and diquark source. We will also 
attempt to isolate the superfluid (and non-superfluid) states of the 2+ld NJL 
model with respect to chemical potential and temperature. The variation of 
superfluidity with respect to temperature will be of most importance. The 
temperature where superfluidity changes from zero to non-zero should be 
the critical point of the system, i.e. the point where vortex and antivortex 
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pairs come together to form the superfluid phase. If this critical point of the 
NJL model coincides with the critical point predicted by the KT theory, this 
would then be strong evidence that the critical phase of the 2+ld NJL model 
could be a gapless thin film BCS superfluid. 



Chapter 9 

Forced Baryon Current Flow 



Chiral symmetry restoration of the 2 + Id NJL model is accompanied by 
the onset of baryon charge density n^, as discussed in Sec. 18.31 in the 
first part of this chapter we will extend this baryon density to a baryon 
3-current where is the timelike component, Jq, of the current. We 
will then implement a twisted diquark source i.e. a periodically varying 
source spanning the whole system. In the final section we will then conduct 
simulations with this twisted source; these will involve comparing the induced 
baryon 3-current with the constant source to one with a twisted source. We 
will show that by introducing a twisted source we introduce a diquark phase 
gradient between the system sites, this then results in a non-zero superfiuid 
current. We will then present the results of the effects on this baryon 3- 
current, as a result of variations in spatial volume, temperature (i.e changes 
in temporal volume) and number of cycles of the diquark source. 
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9.1 Simulating the Baryon 3-Current 

In Sec. 18.31 we saw that the baryon charge density is given by: 

/o = hb =< ^loiJ >= =< ^tr(^-iA') > (9.1) 

This was easily derived, as the NJL lagrangian contains a chemical poten- 
tial term, which allows simple measurements of the baryon charge density. 
This charge density can be denoted as the timelike component of a baryon 3- 
current, which is the conserved Noether current associated with baryon num- 
ber conservation U{1)b, Eqn. ()7.13|1 . The association of a conserved Noether 
current to a global symmetry in continuum field theory has an equivalent 
in Euclidean quantum field theory: known as the Ward-Takahashi identi- 
ties. These identities can be derived for a general action (Using the naive 
fermions): 

•S" ~ ^ ^{i^xlu'ipx+u - i^xlu'ipx-u) (9.2) 
xu 

with a partition function: 



Z = J D^D^e"* (9.3) 
If this partition function is invariant under a simple U{1) transformation: 

V-x ^ i^xc"' i^x ^ ^xe-^" (9.4) 

we can derive the Ward-Takahashi identities, and hence the conserved cur- 
rents of the symmetry. Under the U{1) transformation, the above action 
becomes: 

[e'^^ii'xlui^x+u-i'xlvi^x-u) 
+ e*"(V'x._^7i.?/'^ - ip^+y^yilj^)] (9.5) 
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which can be written as 5 — 5* + 6S, where 

6S f [itpx+ululpx + i'x7u'ipx+u) 

= -iaV~h{x) (9.6) 

where the final line, known as the Ward identity, is written in terms of the 
backwards operator, V~, and Iu{x) is given as: 

^uix) = ^{ipx+ululpx + -ipxlu-ipx+u) (9.7) 

This can be shown to be the conserved current of the U{1) symmetry, by 
writing the transformed partition function as^: 

= J VtP'VtP'il - zaV;/.(x))e-^('^'''^') 
= Z 

=^ < V^h{x) >= (9.8) 

The last line shows a zero expectation value of the Ward identity. This then 

implies that the current /^(x) is conserved. 

Now the conserved current in Eqn. (|9.7|) can be used to determine the 

current components Ji and I2 associated with the U{1)b global symmetry. In 

order to do this we can compare our NJL action, Eqn. (|8.16|) and ()8.19|) . with 

the above action, Eqn. ()9.2jl . so we can write the Ward current components 
^VipVip Vijj'V^' as the Jacobian=l. 
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for our NJL model as: 



1 

lii^) = ::{ViXxXx+i + ViXxXx-i) (9-9) 



where i = 1,2. Now comparing the Iq component with the baryon charge 
density expectation value, Eqn. ()8.4n|) . we can adjust the matrix ()8.38|1 to 
correspond with the Jj components. Then the equivalent matrix for the Jj 
components is given by: 

'~ 'iys'^Kx-i+^^xj ; ^'-^'^ 

where the expectation value of the baryon current density, Jj, are given as^: 
/, =< ij-fiij >=< ^tr{A''X\) > (9.11) 

where i = 1,2. 

9.2 Introducing the Twisted Diquark Source 

Previous simulations, as in [12], were conducted using a constant diquark 
source, 

j^P^P+j^P^P (9.12) 

where j = j = jo- This was sufficient to force diquark condensation and 
derive physically meaningful results in the limit j ^ 0. Such simulations 
have shown, as the Mermin- Wagner theorem states, there is no observable 

diquark condensation for the 2 + Id NJL system. However, the diquark 

^These expectation values are calculated using the Hybrid Molecular Dynamics (HMD) 
code, Appendix El 
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condensate behaves in a non-analytic way, which could imply the existence 
of a type of superfluid phase at high density (discussed in Sec. l8.3.2|) . Such a 
two dimensional superfluid state, without an obvious diquark condensation 
phase, can be associated with the binding [and unbinding] of vortex pairs 
that leads to the critical phase, as in the XY model (Appendix iBj) . For our 
NJL system the vortices can be associated with the phase of the current 
carrying superfluid particles, i.e. the phase of the diquark pairs: 

<i)i) > (f) = ip{x)diquark = (fo exp{i9{x)) (9.13) 

Where 9{x) is the phase corresponding to the external diquark source; where 
for a constant real source, as used in [12], the phase is 6* = 0. 
We can now substitute Eqn. ()9.13|) into the current density: 

Ir^{ipVip*-ip*Vip), (9.14) 

which will give us the baryon current density over the x-y plane, i.e. the 
components Ji and I2: 

1 = TW(f), (9.15) 

where T is defined as the constant of proportionality between the current 
density and the gradient of the phase of the superfluid components, T = 
(0o)^7 and is known as the helicity modulus (Appendix IB|) [38j . From this we 
can see that using the constant source above, or a complex constant diquark 
source: 

j=joexp{ie) j = joexp{-ie), (9.16) 
where ^ is a real constant, we have no gradient over the x-y plane: 

V9 = (9.17) 
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This zero gradient leads to a zero current density. For this reason, we must 
implement a source that generates a non zero gradient over the x-y plane, 
which would then lead to a spatially varying diquark wave function, which 
in turn implies a non-zero current density. If we introduce a "twisted" (or 
periodic) source: 

j = jo exp{ie{x)) j = jo exp{-ie{x)) (9.18) 

where the phase, dependent on spatial position, is given as^: 

e{x) = 271^ + 2tt^ (9.19) 

then, depending on the direction of the twist, this leads to a non-zero current 
density, given by inserting ()9.19|) into ()9.15|) : 

h = T— l2 = T— (9.20) 

So, in an attempt to force a superfluid current flow, we implemented a gen- 
eralised version of EH 



j = ioexp(z27™i^) j = joexp(-z27m2^) (9.21) 

Li 1j2 

where ni and n2 are the number of cycles in the source spanning the lattice. 
Thus, similarly to ()9.2fl|) . the baryon current should become: 

r.^2n ^ ^„27r , , 

h = r—n, I2 = T—n2 (9.22) 

Li L2 

Also, as the baryon charge density would only depend on the magnitude of 
the diquark source, the signal for Jq should be unaffected with this twisting 
of the source. 



^Li and are the x-space and y-space dimensions of the lattice. 
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9.3 Simulations of the 3-current 
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Running simulations, as is often the case, is a time consuming task. Each run 

of the Hybrid Molecular Dynamics (HMD) code'', with particular magnitudes 

of diquark source \j\ = jo and density (chemical potential) /i, generates 

results for one value of /. So, in order to create a plot of baryon current 

as a function of j with 9 points, we had to run 9 simulations - for one 

particular value of chemical potential. Then, to plot the extrapolated value 

of baryon current (as j ^ 0) as a function of /i, with 10 points, would 

mean we need to run 90 simulations. Taking into account the fact that these 

simulations were conducted with variations in spatial dimensions Ls (to study 

finite volume effects) and variations in the temporal dimension Lt (to study 

effects due to temperature variation), we were looking at over 1000 runs. So 

in an attempt to eliminate unnecessary runs we conducted our study in two 

parts. This section represents the preliminary study of the superfluid phase. 

Here we start off first by testing the baryon three current with the twisted 

source. Then we go on to study the general thermodynamic behaviour of 

the superfluid state; that is, we analyse how the baryon current behaves at 

high and low temperatures and chemical potentials, without running detailed 

simulations in the intermediate transition period. In this way, if we see a 

jump in the superfluid current from high temperature to low temperature, 

we could then go and study the transition period, and relate it to the KT 

transition as described by the XY model. Then in a similar way we could 

study the effect of the current with respect to changes in chemical potential. 
''See Appendix El 
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Physical Units 

The simulations in this work were conducted with a strong coupling g"^ = 2.0 
corresponding to Sa = 0.71, as in [12]; where the result Sa is the dimen- 
sionless measure of the constituent quark mass. In the same way our mea- 
surements of the helicity modulus, from the HMD simulations, will be in the 
dimensionless form T = Tga, determined by a straight line extrapolation 
of the baryon current against diquark source. Therefore, in order to obtain 
physically meaningful results we will need to take the ratio: 

i = = 071 i^-^^) 

This will become necessary in Chapter. when we relate our measured 
value of T to temperature. 

9.3.1 Constant and Twisted Diquark Source 

In Fig. 19.11 we have plotted all the components of the baryon 3-current as a 
function of the constant diquark source j = j = jo- In this plot it can be 
seen that the only non zero parts are the baryon charge density Jq. This is as 
we expected, as the constant source produces no gradient, hence there is no 
superfluid baryon current flow, i.e. Ii = I2 = 0. In Fig. 19.21 we have plotted 
the non-zero components of the baryon 3-current with the twisted diquark 
source ()9.21|) . We conducted two sets of simulation, the first was with the 
source only twisted in the x direction, and the second with the source only 
twisted in the y direction. We found, as expected, that for a twist in the 
X direction we had a super current signal Ji 7^ and for a twist in the y 
direction we had a signal I2 ^ 0. These simulations were conducted for the 
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Figure 9.1: Real and Imaginary parts of the baryon 3-current, as a func- 
tion of constant diquark source, j. (Lattice Size=16 x 16 x 64) 



same magnitudes of the diquark source, so that both simulations gave the 
same signal for the baryon charge density Iq = ub- It is also seen that the 
signals generated for the super currents /i and I2 are identical; this is obvious 
from eqn. 19.221 if we have Li = L2 and rii = n2- Therefore in the rest of the 
simulations we only conducted simulations for a twist in the y direction, and 
assumed that a twist in the x direction would give identical results. 

9.3.2 Twisted Source with n-Cycles 

We conducted simulations to test how the baryon current, I2, changed as we 
varied the number of cycles (or twists) in the diquark source, 77-2. This was 
done for high and low chemical potential = 0.8 and 0.2, and the results 
have been plotted in Fig. 19.31 You can note that as you increase the number 
of cycles for the low chemical potential, fi = 0.2, the oc behaviour, 
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Figure 9.2: Non-zero components of baryon 3-current, as a function of 
twisted diquark source j^, with rii = n2 = I. (Lattice Size=16 x 16 x 64) 




Figure 9.3: Baryon current, component I2, as a function of twisted diquark 
source, j2, for various values of = n2- (Lattice Size=16 x 16 x 24) 
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which is more prominent in Figs. 19.61 and 19 .71 becomes more pronounced. 

However, at high /i, the behaviour shifts from an approximate linear fit 
(/^ oc j) to a parabolic fit (/^ oc j^). This shift then leads to a zero extrapola- 
tion of the baryon current as j — 0. So as you increase the number of cycles 
the baryon current ceases to flow. This shift and the resulting destruction 
of the baryon current can be considered due to finite volume effects. From 
Eqn. (I9.22jl it can be seen that increasing n from 1 to 2 cycles has an effect of 
decreasing the wavelength from L to L/2. Therefore, increasing the lattice 
size would counteract the effect of increasing the number of cycles. This can 
be seen in Fig. 19.41 where we have plotted the baryon current, I2, against 
diquark source with n = 2 cycles. We have also plotted the a baryon current 
for n = 1/2. However, due to periodic boundary conditions the results for 
non-integer values of n do not coincide with those of the physically relevant 
current signals, corresponding to integer values of n. This is seen clearly 
from the plots for Lg = 16,n = 0.5 and Lg = 32,n = 1, in Fig. 19.41 

It can also be noted that by using different size spatial dimensions we 
obtain very different baryon currents, so the baryon current is dependent on 
the spatial lattice size. In order to obtain a result that is independent of 
spatial dimensions, we need to work with the helicity modulus, which will 
only vary depending on spatial size as a result of finite volume effects. The 
use of the Helicity modulus will be discussed in further detail in Sec. 19.41 

9.3.3 Effects of Temperature and Density Variation 

In Fig. 19.51 we have plotted the baryon current, /2, as a function of diquark 
source, for a fixed spatial volume (L^ = 16), two values of chemical potential 
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Figure 9.5: Baryon current, I2, as a function of twisted diquark source, j, for 
various time dimensions and fixed space, Lg = 16. 
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and various time dimensions (Lt). At low chemical potential, /i = 0.2, it can 
be seen that there is not much variation in I2 with changes in temperature 
(T = l/Lf). However, at high chemical potential, /i = 0.8, the current can 
be seen to collapse onto a straight line as you decrease temperature. By 
= 128 it becomes clear that I2 extrapolates to a non-zero value as j ^ 0. 
From these results we decided to conduct further simulations using Lt = 64, 
as this gives us results comparable to Lt = 128, however in a third of the 
CPU time. 

In Fig. 19.61 we have plots of various /i at low temperature {Lt = 64). It 
can be seen that the current behaves in a fundamentally different way at 
high density (/i > 0.68) compared to low density (/i < 0.65). As we go from 
high to low /i the negative curvature of the high density phase shifts sharply 
to a strong quadratic type, which can then be extrapolated to I2 = 0. In 
this way we can determine the critical chemical potential for the onset of 
superfluidity as 0.65 < fic < 0.68, which is also the region of chiral symmetry 
restoration i.e. /ic ~ 0.65 (Fig. 18. 2|) . On the other hand, it can be seen from 
Fig. 19.71 that at high temperature {Lt = 4), all the chemical potential curves 
follow a quadratic extrapolation to /2 = at j = 0. So, at high temperature 
no superfluid state exists. Also, even though we have a possible linear 
extrapolation for the low temperature and high density phase, as seen in 
Fig. 19.51 (for Lt = 128), we do not have a systematic method of extrapolating 
to j ^ for /i > /ic in general, therefore we are unable to make any decisive 
statement about the nature of the transition. 
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Figure 9.6: Baryon current, J2, as a function of twisted diquark source, j, for 
various chemical potential (/i) at low temperature (Lf = 64). (Lg = 16) 




Figure 9.7: Baryon current, I2, as a function of twisted diquark source, j, for 
various chemical potential (/i) at high temperature {Lt = 4). (Lg = 16) 
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I2 as a Function of j 

The behaviour of the baryon current, I2, as a function of source, j, can be 
summarised as: 

h{3, /i, T) = /2(0, /i, T) + r)j + T)f (9.24) 

At low temperature, and for high /i (as for the case fi = 0.8 in Figs. 19.51 and 
19. 6 j) we have: 5 = and A 0, and a possible straight line fit extrapolates 
(as j 0) to 12(0, /i) 7^ 0. However, for low /i (as for /i < 0.65 in Fig. 19. 6j) 
we have: A = and B ^ Q which leads to a fit, extrapolating (as j 0) 
to 12(0, /i) = 0. On the other hand, at high temperature this behaviour 
changes (Fig. 19. 7j) . here the extrapolation is of the form for all values of 
/X, i.e. 12(0, /i,T) = y4(yU, T) = 0. The extrapolation shift from j to can 
be understood if we have different wave functions for the current carrying 
superfluid particles at different chemical potentials: 

0(f) = (?!)ojoeXp(z6'(x)) yU < /ic 

0(x) = 00 exp(z6'(x)) /i > /ic (9.25) 

The linear jo dependence, at low chemical potential, is expected due to the 
fact that there should be no diquark condensate for /i < /ic; so </> ^ as 
i 0. In this way, using Eqn. ()9.14|1 . the baryon current is: 

I ~ 0oJoV^(f) ^^< 

I ~ <PlVe{x) fi>fic (9.26) 
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and in the limit j 0, we have: 

/ ~ fi < fic 

I ~ (l)l'V6{x) = constant > /j-c (9.27) 

which is the behaviour of the baryon current as seen throughout our simula- 
tions. 

9.4 Finite Volume Effects 

We have, so far, been working with the baryon current, I2, for example in 
Fig. 19.61 we have plotted the current against diquark source. However, this 
is not the most appropriate variable to work with, due to its dependence on 
spatial volume, as we saw briefly in Sec. 19.3.21 This dependence on spatial 
volume can be seen clearly when we conduct simulations for a flxed tem- 
perature, but differing spatial sizes: in Fig. 19.81 it can be seen that, at low 
temperature, the baryon current differs hugely for the three different spatial 
volumes Lg = 16, 24 and 32, and any linear extrapolation of the baryon 
current with j ^ will give very different results. However, at low jj, and/or 
high temperature, when the extrapolation is of the form, finite volume 
effects do not matter too much as all extrapolations tend to J2 — > 0, Figs. 19.81 
and EH 

9.4.1 Working with the Helicity Modulus 

We can overcome effects of spatial dependence by working with the helicity 
modulus, T, which, as seen earlier, is the constant of proportionality between 
the baryon current and the gradient of the source's phase, Eqn. ()9.15|1 . In this 
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Figure 9.9: Baryon current, I2, as a function of twisted diquark source, j, for 
various space dimensions, /i = 0.8, and high temperature, Lf = 12. 
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Figure 9.10: Helicity modulus, T, as a function of twisted diquark source, j, 
for various spatial dimensions, /i = 0.8 and low temperature, Lt = 64. 

way we are able to represent, in a volume independent manner, the superfluid 
behaviour of the system. So, when we plot T = l2Ls/2n against j for varying 
spatial volumes, L^, but constant temperature and chemical potential, any 
variations in the plots could be put down to finite volume effects. 

Low Temperature Finite Volume Effects 

In Fig. I9.1()l we have plotted T against j for low temperature, Lt = 64; this 
is only done for high yU, as for lower values of /i the extrapolation is always 
of the type, which implies the current, and thus, the helicity modulus 
extrapolates to zero as j 0. We can note that the dip in T, at low j, 
decreases as the spatial volume increases; thus, if this dip is a finite volume 
effect it should tend to a straight line in the infinite volume. In order to 
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Figure 9.11: Helicity modulus, T, as a function of 1/Ls, for various j, /i = 0.8 
and low temperature, Lt = 64. 



analyse this dip at larger volumes, we extrapolated T to the limit of Ls —>■ oo 
by plotting T for each value of j against 1/Ls, Fig. l9.1l( and then the infinite 
volume value of T was deduced by extrapolating the plots l/Ls — > 0. The 
infinite volume plot of T is also plotted in Fig. 19.101 and it can be seen 
that the dip vanishes in this limit, hence this effect can be put down to 
finite volume effects. From the low temperature, infinite volume plot, we 
have conducted an extrapolation of T for j ^ 0, which should represent the 
dimensionless value of T at low temperature: T = 0.1413 ± 0.0014. 

High Temperature Finite Volume Effects 

Using the same techniques as before we studied the superfluid behaviour at 
high temperature, Lt = 4. In Fig. 19.121 we have plotted T for /i = 0.8, = 4 
and, again, for various spatial volumes. As before we extrapolate to the 
infinite volume limit, where the extrapolation lines are shown in Fig. 19.131 
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Figure 9.12: Helicity modulus, T, as a function of twisted diquark source, j, 
for various spatial dimensions, /i = 0.8 and high temperature, Lt = 4. 
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Figure 9.13: Helicity modulus, T, as a function of l/L^, for various j, 11 = 0.8 
and high temperature, Lt = 4. 
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In Fig. l9.12l we have also plotted the curve for Lg oo. From this we see that 
at high temperature and in the infinite volume limit the T oc behaviour 
does not change significantly, if anything it becomes more pronounced at 
lower values of j. From this extrapolation we find T = 0.00012 ±0.0012 ~ 0. 

Behaviour at Intermediate Temperature 

We have shown that at low temperature there are significant, but control- 
lable, finite volume effects, which are removed by extrapolating to the infinite 
volume limit. In the infinite volume limit the low temperature linear be- 
haviour, T(j) oc T(0) + Aj, becomes obvious and straight line extrapolations 
for j — 5> lead to T 7^ 0. At high temperature, where we expect a T(j) oc 
behaviour, working on a finite volume does not have much of an eff'ect, that 
is the finite volume plots do not differ from the infinite volume extrapolation. 
So, we can state that the superfluid current drops from a non-zero value at 
low T to a zero value at high T. 

However, if we attempt to analyse the intermediate temperature region 
we encounter some difficulties. In Fig. 19.141 we have plotted T for /i = 0.8, 
Lt = 12. The extrapolation of this, to the thermodynamic limit, is given 
in Fig. 19.151 We see that the various volume plots have an "s" shaped 
form, which remains even in the infinite volume limit. In this case we could 
implement a extrapolation for the first four points, however this does not 
seem very conclusive as the last 4 points lie well outside any behaviour. 
On the other hand, in Fig. 19.161 where we have again plotted a curve for 
Lt = 12, but now we have used a larger volume, (L^ = 40), and obtained 
further points at low j, a extrapolation does seem appropriate. From this 
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Figure 9.15: Helicity modulus, T, as a function of l/L^, for various j, fj, = 0.8 
and high temperature, Lt = 12. 
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Figure 9.16: Helicity modulus, T, as a function of twisted diquark source, j, 
for Ls = 40, /i = 0.8 and high temperature, Lt = 12. 



extrapolation T can be deduced to be T = 0.00016 ± 0.0014 ~ 0. 

A similar situation occurs as you decrease the temperature further. In 
Fig. 19.171 we have a plot of Lt = 24, and its extrapolation lines in Fig. 19.181 
At this temperature we can see that it is not obvious if this plot falls under 
a linear extrapolation (giving T 7^ 0) or a parabolic extrapolation (giving 
T = 0). At low j the curve must tend to T = hence it must follow a 

curve to the origin. But, more plausibly, we could apply a straight line 
extrapolation to the last 5 points. Even in the infinite volume limit it does 
not obviously fall under either extrapolation criteria. 

In the next chapter we will go on to investigate the intermediate tempera- 
ture range. We will attempt to deduce some criterion for the extrapolation of 
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Figure 9.18: Helicity modulus, T, as a function of l/L^, for various j, fj, = 0.8 
and temperature, = 24. 
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these temperatures, which will then help in isolating the critical temperature 
i.e. the temperature at which T falls from the non-zero superfluid value at 
low T to a zero value in the non-superfluid phase at high T. 



Chapter 10 

Thin Film Superfluid Dynamics 



In the previous chapter it was shown that the baryon current-density has a 
strong transition in the region 0.65 < fic < 0.68 - the region corresponding 
to chiral symmetry restoration. Also, working with the spatially indepen- 
dent helicity modulus, T, it was shown that at low temperature, Lt = 64, 
and high chemical potential, = 0.8, the helicity modulus extrapolates to 
T = 0.1413±0.0014; and at high temperature, Lt = 4, it extrapolates to zero. 
However, at intermediate temperatures, for example Lt = 12 and Lt = 24, 
extrapolation techniques are not obvious - even in the infinite volume limit. 
In this chapter we will investigate the behaviour of T in this intermediate 
region, and attempt to provide some extrapolation criteria for this puzzling 
domain. The aim of this chapter will be to try and isolate the critical tem- 
perature and see if it corresponds with the predicted T^t temperature, given 
by the theory of Kosterlitz and Thouless (KT). 
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10.1 Tkt given by Kosterlitz and Thouless' The- 
ory 

If the phase transition, for the 2+ld NJL model at high density, is described 
by the KT theory (Appendix [HI) then the critical temperature should be 
given by Eqn. |R23| . or in lattice units: 

12 1 ^ 12 1 , , 

aLt = —- = -— Lt = -—- = --, (10.1) 

where T = Tga is the dimensionless measured value. Using our value of T 
we find: 

2 1 

4.51 ±0.44 (10.2) 



vr 0.14129 ±0.00138 
Therefore, according to the KT theory the critical Lt should be approximately 
6, which is the nearest even lattice length where the fiuid phase jumps to that 
of a superfiuid state. This implies that everything above (and including) this 
critical length (temperature) should be in the superfiuid phase (the bound 
vortex phase)i.e. T 7^ 0. This suggests that the Lt = 6 data should follow a 
linear extrapolation to the origin. However, from our previous findings it was 
shown that Lt = 12 follows more of a extrapolation than a linear one, but 
even this is not very clear. We will now investigate the intermediate domain, 
ranging from Lt = 6 — 64, in order to determine, whether or not, the KT 
critical point exists in our 2 + ld NJL model. 

10.2 Extrapolation Criteria and the Critical Point 

In order to investigate the change in T from the high temperature (T(j) oc j^) 
phase to the low temperature linear (T(j) oc T(0) ± Aj) phase we have 
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Figure 10.1: T as a function of j for various Lt (/i = 0.8 and Lg = 32). 

plotted, in Fig. 110. IL T against j for temperatures ranging from Lt = 2 to 
Lt = 64. In this plot we can see a steady transition from the obvious 
form at Lj = 2 to an "s" shaped form between Lj = 6 — 24 and a straight 
line form above Lt = 24 (which dips to the origin at low j). From this 
we cannot really distinguish any real, abrupt, phase change. So, to obtain 
a different perspective we have plotted the same data as T against Lt for 
the various j values in Fig. 110. 2[ and in Fig. 110.31 we have a more detailed 
plot of the j = 0.025 data. From Fig. 110.21 we can see that T decreases 
steadily from Lt = 64 to Lt = 24, for j > 0.025. However, below Lt = 24 it 
drops drastically to T = 0. For j = 0.025, in Fig. 110.31 it can be seen that 
T follows almost a straight line. In Fig. 110.41 we have conducted a straight 
line extrapolation to T = 0, from which we found Lt = 5.54 ± 1.66. Also, in 




Figure 10.3: A detailed plot of T as a function of Lt for j = 0.025 (/i = 0.8 
and Lg = 32). 
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Figure 10.4: T as a function of Lt for j = 0.025 with extrapolation line 
y = aox + ai (ao = 0.00192 ± 0.0005 and ai = -0.01062 ± 0.0015); (/i = 0.8 
and Lg = 32). 



Fig. ll0.5l we have conducted another extrapolation, of the same data, but now 
with only the first 18 points; from which we find: Lt = 7.24±0.9. From these 
two fits we have the average: Lt = 6.4 ± 1.9. This value of Lt is comparable 
to the one derived in Eqn. ()10.2|) . Therefore, we could conjecture that the 
critical KT temperature could be described by the point where T intersects 
the Lt axis for a fixed small j. In order to confirm this we would need to 
conduct further simulations using smaller j (e.g. j = 0.01) and maybe larger 
volumes (in order to reduce finite volume effects). 

10.2.1 Finite Volume Effects at Low Temperature 

Even though we have located a possible candidate for the critical temper- 
ature, it still does not provide any obvious extrapolation criterion for the 
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Figure 10.5: T as a function of Lf for j = 0.025 with extrapolation line 
y = aox + ai using only 18 points (oq = 0.002124 ± 0.00025 and ai = 
-0.015367 ± 0.0009); (/i = 0.8 and L, = 32). 



range = 6 — 24. 

Our next task will be to find a definite lower limit for the superfiuid state, 
i.e. the smallest Lt point for which infinite volume extrapolations eliminate 
the dip in T at low j. As for Lt = 64, Fig. I9.1()[ we conducted similar 
extrapolations for T down to Lt = 48. It was found that the dip at low j 
disappears for extrapolations down to Lt = 56, Fig. 110. 6j these simulations 
were only done for 4 points, at low values of j, in order to reduce time. 
So, down to and including Lt = 56 the dip at low j vanishes in the infinite 
volume limit. However, below this point the dip becomes more and more 
pronounced, even in the infinite volume limit. Fig. 110.71 As can be seen for 
Lt = 24, Fig. 19.171 the dip eventually becomes something that you cannot 
just ignore. 
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Figure 10.8: T as a function of l/L^ for = 56 (/i = 0.8 and = 32). 



As in the previous chapter, the infinite volume limit extrapolations for 
Lt = 48 and Lt = 56, Figs. 110.61 and 110. 7L were done by a straight line fit of 
T against 1/L,, Figs. lUTSl and IuTqI 

10.2.2 Critical Temperature Range 

In a conservative conclusion, the only phase range we can be sure about 
is that above Lt = 56 and the one below Lt = 4. For Lt > 56 infinite 
volume extrapolations lead to a perfect straight line, and a non-zero value 
of T, implying superfluidity. For Lt < 4, even without extrapolation, the 
curves are parabolic, Fig llO.H thus, they could be associated with the non- 
superfluid state (T = 0). This then leaves us with an uncertain range from 
Lt = 6 — 54, which is huge! However, as I stated before this is a conservative 
conclusion. As we deduced through the plot of T against Lt for j = 0.025 
(Fig. [TIE3I) the critical point could plausibly occur at Lt ~ 6, in which case 
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Figure 10.9: T as a function oil/Ls for Lj = 48 (/i = 0.8 and Ls = 32). 



all extrapolations above this point would be of the type T(j) oc T(0) + Aj. 
This could very well be the case for 1/T above Lt = 24, even though we can't 
explain the dips at low j. This then leaves us with the range = 6 — 22 as 
the uncertain domain, due to the fact that these curves follow a "s" shape 
which could be either a linear or parabolic extrapolation. 



Chapter 11 
Summary 



In Chapter ini the calculation of baryon charge density (ub) was extended to 
the conserved baryon 3-current where Ub was the timelike component, 
Jo, of the current. Then, to force a current flow we implemented a twisted 
diquark source i.e. a periodically varying source, spanning the whole system. 

We then investigated the current's behaviour with a change of density 
/i, for two temperatures^: ^ = 0.565(1) and ^ = 9.043(1), correspond- 
ing to lattice lengths = 4 and Lt = 64 respectively. We found that a 
superfluid phase does exist at high density (/i = 0.8) and low temperature 
^ = 9.043(1); but for high temperature, ^ = 0.565(1), the current is always 
zero. Figs. 19. 6| and 19.71 For the low temperature case we found that the cur- 
rent density has a phase transition in the region 0.65 < fic < 0.68, Fig. 19.61 
This transition region corresponds to that of chiral symmetry restoration in 

These findings show that the low temperature, high density, phase of 

the 2+ld NJL model is a thin film superfluid, and the transition from this 

superfluid phase at high density to a non-superfluid phase at low density is 
1 Given as the ratio LjT = ^ = ^. 

la 1 
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of 1st order. 

Initially results were presented in a spatially dependent form, so accord- 
ing to the size of lattice used different results were obtained for the current. 
We later showed that it is better to work with the spatially independent he- 
licity modulus, T, defined as the constant of proportionality I2 = T2ti / L2. 
The non-zero value of the helicity modulus was then determined in the ther- 
modynamic limit for low temperature = 9.043(1)) and high chemical 
potential (/i = 0.8) as T = 0.1413(14) . We also showed that this falls to 
~ at higher temperatures (^ = 0.565(1)). Therefore, the ratio T/S that 
represents the superfluid phase of the 2 + Id NJL model at high density is 
given as = 0.199(1) ~ 0.2 

Once we had found the presence of a superfluid phase, next, we wanted to 
test whether this superfluid behaviour was due to the binding and unbinding 
of vortices, as described by the KT theory for 2-dimensional critical systems. 
In Chapter [101 we started by calculating the critical temperature predicted 
by the KT theory. Using our value of T this critical temperature was found 
to be ^ = 0.637 ~ |, which implies a lattice length Lt = 4.51(44). 

We then showed that by plotting T against Lt for low j {j = 0.025 in our 
case), we obtained a very good straight line flt that intersects the Lt axis at 
a point near the KT critical length prediction, Fig. ll().4l and llO.,^ Through a 
straight line extrapolation (T oc ALt + B) this critical point was found to be 
at Lt ~ 6(2). This point corresponds to the KT prediction within error. Also, 
this critical transition point seems to be of 2nd order. However, in order to 
make any conclusive statements about the critical point, further study is still 
needed, maybe using larger volumes, Lg > 32, and smaller diquark source 
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values, j < 0.025. 

We also showed that at low temperature, ^ > 7.9 (or in terms of lattice 
size Lt > 56), the dip of T at low j can be put down to finite volume effects, 
which disappears when we extrapolate to the infinite volume limit. However, 
even though the plots between Lt = 24 — 54 do follow straight lines for high 
j, similar infinite volume extrapolations do not eliminate the dips at low j. 
The domain between Lt = 6 — 22 is even more puzzling. In this area we 
cannot easily apply a straight line extrapolation, as the curves of T against 
j are parabolic at low j. 

So, in order to construct a conclusive phase diagram for T against Lt, we 
not only need to confirm the location of the critical point, but we also need 
to investigate and explain the strange behaviour of T in the intermediate 
temperature range. 



Appendix A 

Schwarzschild Spacetime 



The 'Schwarzschild geometry' describes a spherically symmetric spacetime 
outside a star, and its properties are determined by one parameter, the mass 
M. The Schwarzschild metric, in spherical polar coordinates, takes the form: 



1 F 
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^2 





V 








— sin(6') / 



(A.l) 



This type of spacetime geometry is said to be static, due to the fact that (i) all 
metric components are independent of t, and (ii) the geometry is unchanged 
by time reversal, t —t^. 

A.l Geodesic Equations 

The equations of motion in Schwarzschild spacetime can be derived by using 
the covariant form of Newton's second law of motion. 

Dp'' 



Dt 
dp'' 

dr 



~^ _|_ TiP 

1 ' 1, 



up 



dr 



(A.2) 



space time with property (i) but not necessarily (ii) is said to be stationary i.e. a 
rotating star/black hole 
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where we have p^^ = and for free fall we have F^' = 0, then the 



equations of motion are: 



(fx^ ^„ dx" dxP , . , 



Now, using this and the appropriate Christoffels we are able to determine the 
equations of motion. Alternatively, by using the Schwarzschild Lagrangian 

L = - —)di - (1 - — )-Mr' - r'^di - (r2sin2 0)d'/] (A.4) 



2' 



we find 



dL dL , 2M,. 

d dL , 2M,.. 2M . 

-AX^ = ^ )t + —rt 

dX at r 



dL M ^ M\ 2M 2.2 n2 V2 . 2n 

— — = —t H (1 ) r — rd —rep sm 

dL 2M 

-wr- = -1 r 

or r 



d dL , 2M,_T 1M , 2M. 

-jTjrr = -1 r+— 1 ^ 

dX or J. J.Z y. 



2 . /) d^ 2/j 

— — = —r smticosUd) — - = —r 6 

de ^ d9 
d dL o-A 

dX de 



dL n dL 2 • 2nl2 

—— = — - = —r sm 6(p 



— = —2rr sin^ Od) — 2r^ sin 9 cos Od) — sin^ i 

dXd6 



(A.5) 



(A.6) 



(A.7) 



(A.8) 



APPENDIX A. SCHWARZSCHILD SPACETIME 147 



and substituting into the Euler equation 

dL d dL 



we have the equations of motion for Schwarzschild spacetime 



(A.9) 



2M,.. 2M . , 
= (1 - — )t+ — ft = (A.IO) 



= — - - sin^ ^ + (1 

= -r2sin^cos^</)^ + 2rf^ + r26' (A.12) 
= 2rf sin^ + 2r^ sin 6^ cos 6^0 + sin^ 6*0 (A. 13) 



A. 2 Riemann Components in the Orthonormal 
Frame 



The spacetime line interval, Eqn. p.lljl . can be rewritten using tetrad trans- 
formations (discussed in the next section): 

OAT OA/ 

u' = {l-—)^dt uj' = {l-—yUr J = rde uj'^ = r sm{e)d(t) 

(A.14) 



as 



ds^ = - icu'-? - icu'? - (cu^? 



(A.15) 



This now has the Minkowski metric 



/ 1 \ 
0-100 

0-10 

\ -1 / 



(A.16) 
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Noting that the metric can be written as: 

dgf,u = ujf,u + uj,yf, (A. 17) 

and using the fact that: 

dg,. = ^dx- = (A.18) 

we then have|2]: 

ujf,u = -^^u^l (A. 19) 

which impUes u is antisymmetric, i.e. it has only six unique components and 
^^fj.u = for /i = z/; and we have the conditions: 

u;^ = uj^, u} = -ui (A.20) 

We can now write the exterior derivatives^ of 



1, 2M,_i,2M. 
M ,^ 2M. 1 



du^ = -(1 —)-2(—)drAdt 



r 



2 



1 )—^Lu''Au' (A.21) 



du^ = (A.22) 

duj'^ = dr Ad9 = ^{l-'^)^u;'^ Au^ (A.23) 

duj'^ = sm{9)dr A d(p + r cos{(j))d9 A dcj) 

= l(l-^)VAc.^ + ^o;^Aa;^ (A.24) 



Now, using Cartan's equation 

du^" = u'' Au^^ + n" (A.25) 

^The exterior derivative: d = dx^^^ acting on a 1-form A'^ = A{x)dx^ gives dA'' = 
dx^' A dx"^^, where dx^ A dx'' = for n = 
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with zero torsion (^2^ = 0), we can write 

du^ = uj"^ A iul + iu^ A Lul + Lu''' A (A.26) 

dcu" = Aul + uj^ Aul + u^ Au'^ (A.27) 

du'^ = Auf + uj'' Au'l + uj^ Aul (A.28) 

da;* = uj' A ujf + uj'' A ujf + 00^ A uj^ (A.29) 



Comparing Eqns. (TOTTl - irOl with Eqns. (r06ll - (r09l) we find the six 
unique components of cuf^ as 



, M, 2M, 1 , M , , . . ^ 

cc', = ^(1 )~-2uj' = —dt = ul A.30 

. 1, 2M,i . , 2M,i ,^ , , . X 

CJ,^ = -1 2,^^= 1 Y^^Q=_^^r 

ujf = -{l-—)-2uj'f' = sm{e){l-—)-2d<P=-u;l (A.32) 

ry ry 'y ^ \ / 



cot(^)^^ = cos(^)ci0 = -o;J (A.33) 



r 



oj^'' = u)l = Q (A.34) 
a;f = o^i = (A.35) 



Now, using: 

i?^; = dojf; + < A < (A.36) 

we can write the six independent components of the Riemann tensor: 

Rl = dujl + cj^ A cj" = dujl + ujIa uj^ + iulAujf 
2M . . 2M 



dr Adt = — ^o;"" A cc;* 
2M 



^3 



^rrt Rtrrt 



(A.37) 
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2M, i2M , ,^ M , e 



— Rf)r; 



M 



■Orre 



(A.38) 



Rf = dujf + iotAco^ = dujf + ufAiul + iu^Aiu^. 

2M 1 M 2M 1 
= cos(6')(l )2d6' A rf0 + sin(6') — (1 r^dr A 

2M 

+ cosmri f2d(t)Ade 

r 

cot(e), 2M,i e <6 M , , 



M 

3 



cot(^) 2M 

r 



r 



(1 —)^^^ Auj= —OJ' A UJ^ 

M 



(A.39) 



Ri 



dut + ujtAu;^ = dut + uf Auji + A 



iUn 



-sm{9)d9 Ad^-sm{9){l 



2M 



)dcf> A dd 



J Auj^ + )J Au^ = -u^ A uj^ 



ry'2i iy'2i ly 



Ri 



-R. 



2M 



(A.40) 



R' 



dujl + uj1a lu^ = dul + ujI A ujI + ujI A ul 



M , 

- 

^tet - 



2M. 1 



)^d9 A dt 



M 



-co" A cu' 



r 

—Retet 



M 



(A.41) 
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Rt 



dujf + uJi^ uj^ = duf + iuf A iul + iu^ A uol 
^(1 - — )^ smie)de Adt = ^uj'f' A uj' 



t(j)t 



-R 



M 



(A.42) 



A. 3 Tetrad Transformation 



In Eqn. (|A.14|) we transformed from the coordinate frame to an orthonormal 
frame. This is achieved by using the tetrad of the form: 



/ 

F-5 

r 

V rsin(^) / 



(A.43) 



This takes a 4-vector, given in the coordinate frame, and maps it to the 
equivalent in the orthonormal frame. We can then determine the inverse 
tetrad, i.e. a tetrad that takes vectors in the orthonormal frame and maps 
them to the coordinate frame. We find this by using the inverse of (|A.43j) . 



la 












i 

r 











(A.44) 



which has the property: 



St 



(A.45) 



Therefore, (jA.44|l takes vectors in the orthonormal frame and maps them to 
the coordinate frame. 
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A. 4 Riemann Components in the Coordinate 
Frame 

In Sec. IA.2I the unique components of the Riemann tensor were calculated 
in the orthonormal frame. In order to determine the components in the 
coordinate basis we will use the transformation tetrad given in Eqn. r04l 
Before we can us the tetrad, (lA.44j) . we need to adjust it so it is able to take 
covectors in the orthonormal frame, rather than vectors: 

/ F3 \ 





r 
V Or/ 



(A.46) 



which is just Eqn. ()A.43|1 . Then using this we have: 

^ (coordinate frame) = (e~^)|JAa(orthonormal frame) 

Now using this tetrad we can transform the orthonormal Riemann tensor 
components as: 



RfMuap = (e ^)1^{e ^)|^(e ^)^(e ^YpRabcd 
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Now using the components given in Sec. IA.2I we have: 

Krro = {e-')l{e-'yr{e-'):{e-')lRerre = rF-'^FM-^) = 
Krr, = (e-^)J(e-^);(e-i);(e-i)X-0 = ^F-^fM'^) = 

Km = {e-Xie-yAe-reie-X,Rem = rF'^rFH-^) = -^ 
Kt,t = ie-X{e-yAe-X(^-XtKt,t = rF'^rFk-^) = -^ 

(A.47) 

Therefore, the Riemann tensor components in the coordinate frame are given 
as: 

_ _2M MF-' MF-' 

^trrt — ^3 ^err9 " ^ ^<l,rr<t> " ^ 

_ p/ _ ^F p/ _ ^F 



(A.48) 



Appendix B 

XY Model and the KT Transition 

B.l Spin Systems 
B.1.1 The Ising Model 

The most basic spin system is the Ising model, which is generally used to 
describe magnetisation. It consists of spin sites, each taking a value of 
+ 1 or —1 (i.e. up and down respectively). In its simplest form it describes 
a 1-dimensional lattice, usually with the N sites in a chain formation and 
periodic boundaries. Other simplifications include a uniform external field at 
each site. Hi = if, so that the only interaction between spins is that between 
the nearest neighbours, symbolised by J. This system is then described by 
the Hamiltonian: 

N 

-n = Y,HiSi + jJ2 S^Sj (B.l) 

1=1 <hj> 

The partition function of this model is then written as: 

Z = J2^M-m{S,)) (B.2) 

[Si] 

where P = l/^ksT), and the the sum is taken over all possible states i.e. a 
total of 2^ possible spin states. Then, to extract thermodynamic information 
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from the partition function we define the free energy of the Ising model as, 

F = -ksT In Z 

From this we can determine the magnetisation of the system^ 

i=l i=l 

where < Si > is the average spin of the i*'* site. The correlation function, a 
measure of the influence exerted by a given spin Si on other spin sites 5*^ , is 
given as: 

Gij =<i SiSj > — <iSi^<iSj^] 

Since the interaction between spins favours alignment, a nearby spin Sj will 
tend to assume the same orientation as Si; however thermal fluctuations 
counteract this tendency and exert a de-correlating efl'ect. Thus we expect 
some correlation that weakens as the distance between Si and Sj increases, 
also at a fixed distance apart the correlation will be stronger when the tem- 
perature is lower, see Fig. IB.ll [27| . So, when the temperature is low we can 
assume that the majority of spin states are correlated, which in the absence 
of an external field leads to the phenomenon of spontaneous magnetization, 
Mg. For the 1-dimensional Ising model this does not occur for T 7^ 0, but 
there is a trivial phase change at T = 0, i.e. when the thermal energy is zero 
all the spins are aligned. However, in 2-dimensions there is a second order 

phase transition at the Curie temperature T = [23 [211, Fig. IB. 21 

iPor Hi=H we have < S >= < >= |f 
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Figure B.l: Ratio of the correlation function G 
ferent temperatures T (not to scale). 



G(r,j)/G(0), for two dif- 



M 




Tc T 



Figure B.2: Spontaneous magnetization in 2-dimensions against temperature. 
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B.1.2 The XY Model 

The restriction of the Ising model is that the spin vector can only assume one 
of 2 discrete values. Even though for some cases this is fine, a more realistic 
model is the XY model. In this case spin is a continuous vector in a plane 
(for 2 dimensions). This model is characterized by the Hamiltonian: 

N 

-n = Y, H^Sl + JY,S^■S, (B.3) 

i=l <i,j> 

where the spins are classical vectors of unit length aligned in the plane. For 
zero external field we can write: 

<i,j> <hj> 

where $ is the angle between the spin direction and an arbitrarily chosen 
axis, and is a continuous variable. 

In this 2-dimensional model, spin excitations (spin waves) are easier to 
excite thermally than for the Ising model, they are so strong that they can 
destroy long-range order at any finite T; so, in the case of 2-dimensions, 
symmetry breaking for T > does not occur (Mermin- Wagner theorem [29]). 
On the other hand, in 1973, Kosterlitz and Thouless showed that a kind of 
phase transition can exist for the 2d XY model, or more accurately there is 
a critical point at non-zero temperature. It arises not from the long-range 
ordering like spontaneous magnetization but from the excitations of vortex- 
antivortex pairs |3T]|32]. 
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B.2 Phase Transitions in the XY Model 

B.2.1 The Critical Phase in the XY Model 
Vortices 

If we consider small spin variations, $i — $j ^ 1, then we can expand [EH as: 
Hc^Eo + ^jY: (^^ - = ^ E l^^(Or (B.5) 

<i,j> i 

where a is the lattice spacing. This can now be taken to the continuum limit, 
so that $(r) becomes a scalar field, and the sum becomes an integral over 
the 2d area: 

'H-Eo = ^ I dV|V$(r)|^ (B.6) 



2 „ 

In this 2d state of continuous spin states a possible excitation above the 
(perfectly ordered) ground state is a vortex^. Fig. IB. 31 Any line integral 
along a closed path around the center of the vortex will give: 



where q is the vorticity and is the angle to r.^ Now using (|B.6j) and ((BT 
we can find the energy of this vortex above the ground state energy: 

E= r 27vrdr-{^ f = 7rJ\n{-)q^ (B.8) 
Ja 2 r a 

here, R is the size of the system [2ni[2Z|, and a is the lower cut off (as the 
theory diverges as r ^ 0). Now it can be seen that the energy of the vortex 

increases logarithmically with the size of the system. 

^Spin waves are another form of excitation, and are responsible for destroying long- 
range order at any finite T. 
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Figure B.3: Vortex with vorticity |g| = 1 



Thermodynamics of Single Vortices 

A vortex can be put onto any lattice site, thus the number of configurations 
is {R/aY, and then we can write the entropy of the vortex as 

S = 2kBH-) (B.9) 
a 

Therefore, the free energy, F = E — TS, is positive and large at low tem- 
peratures. Thus, it is not very easy to thermally excite single vortices at low 
temperature. But, as temperature increases, the free energy decreases, and 
at a critical temperature Tkt the free energy vanishes: 

kBTKT = ^ (B.IO) 

Then, above Tkt a large number of vortices are thermally excited. 
Vortex- Antivortex "Dipole" Pairs 

Even though it is difficult to excite single vortices at low temperature, it 
does not mean none exist. In the low temperature phase it becomes more 
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favorable to excite pairs of vortices, i.e. a vortex-antivortex pair. To see this 
we start by writing the energy of a system of two vortices, with vorticity qi 
and 

$ = $i(fi)$2(r2) 

R 

E{qi,q2) = -7rJ(gi + g2)^ln(— ) 

a 

- 27rJgig2ln(^^^^^) (B.ll) 

where for a vortex and antivortex pair qi = —q2, thus we have: 

E{q,,q2) = 27rJg2ln(tl^) (B.12) 

a 

From this you can see that the self-energy of a vortex-antivortex pair is small 
when the vortices are close together. Thus, as the energy of a "dipole" pair 
is smaller than a single vortex, it is easier to excite a pair of opposite vortices 
than a single vortex. 

Critical Phase of the XY Model 

We can now see that at low temperature, < T < Tkt, there will be a finite 
density of low energy vortex-antivortex pairs, of zero total vorticity. Then, at 
high temperatures, T > Tkt, there will be a large concentration of unbound 
vortices. 

In this picture the critical temperature is then associated with the un- 
binding of vortex-antivortex pairs; the pairs are bound strongly to each other 
at low temperature, and as the temperature increases, the number of pairs in- 
creases. Then, at the critical temperature Tkt, vortex-antivortex pairs start 
to separate and become free; this is called the KT (Kosterlitz and Thouless) 
phase transition. 
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Correlation Length of Vortex Pairs 

Due to the KT transition the correlation function of the XY model changes 
behaviour as it crosses the critical point Tkt- Above the critical point, when 
the vortices are unbound and free to interact with each other, the correlation 
function behaves in a classic way: 

<SiSj>=e~l, (B.13) 

with the correlation length: 

^ = ^1 

However, below the critical temperature the correlation function shifts to a 
power-law: 

< S,Sj >= (^)t (B.15) 

and the critical exponent becomes: 

1 T 



(B.16) 



e 27rJ 

This can be understood in the following way: at low temperature vortices 
bind into dipole pairs, whose influence on the system is confined to small 
distances, thus the only important large distance interactions are spin waves. 
This phase is then described by the power-law behaviour of the correlation 
function. As the temperature rises the size of the dipole pairs increases, and 
diverges at T = T^t when free vortices appear. These free vortices then 
interact, and in the process, the low temperature dominance of spin waves is 
destroyed. This phase is then described by an exponential behaviour of the 
correlation function. 
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B.3 The KT Transition and Superfluidity in 2d 

One application of the KT transition, i.e the binding and unbinding of vor- 
tices at the critical temperature T^t, is in understanding superfluidity in 
thin fllms. For d > 2 Bose-Einstein condensation occurs in a Bose gas at low 
temperature, and this condensed phase is deflned by the breaking of a U{1) 
global gauge symmetry. However, for systems with d < 2 the breaking of a 
f/(l) global gauge symmetry is prohibited by the Mermin- Wagner theorem, 
i.e. no phase transition to a state with long-range order can occur at finite 
temperatures. On the other hand, topological phase transitions are possible; 
and the binding and unbinding of vortices in two dimensions leads to a type 
of phase transition. 

B.3.1 Superfluids and the XY model 

Our previous calculations, relating to vortices in the XY model can be di- 
rectly applied to the study of superfluids. For single vortices, the XY spin 
angle $ (jB.5j) can be related to the phase of a particle's wave function in the 
superfluid condensate state: 

^ = v/^exp(-i$(t,f)) (B.17) 

where Uq is the number of condensate particles (superfluid density). Then 
the current density of the superfluid is given by: 

Icond = 

2m 

= — noV$ 

m 

= TV$ ^ T = —no (B.18) 

m 
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where m is the mass of the superfluid particle; and T is the constant of pro- 
portionality between the current density and the phase angle, and is known 
as the helicity modulus (T). The phase angle $ can be associated with the 
spin angle as in the XY-model or with a particles phase as in the superfluid 
condensate. 

If we divide both sides of ()B.18jl by the density no, we have the velocity 
of the condensate particles: 

Vcondi^ = -V$ (B.19) 

m 

This can now be inserted into the energy of a superfluid velocity fluctuation, 

n = ^mno j d?r\vcond{^? (B.20) 
where mriQ is the mass density of the superfluid. This then gives us: 

= ^ [ rfV(V$)' ^ T = -no (B.21) 
2 / m 



It can now be seen that this is similar to (lB.6j) if, as stated before, we identify 
the phase of the superfluid condensate with the angle of spin in the XY 
model and hT = J . Therefore, due to this link between the XY model 
and superfluids, the transition from a superfluid state to a normal state for 
a 2d superfluid can be associated with the binding and unbinding of vortex 
pair excitations with opposite circulation. Evidence of this KT behaviour 
is seen in superfluid Helium 4. In 3d there is a continuous decrease of the 
superfluid density with temperature, but in the 2 dimension case, for example 
in thin fllms of ^He, a jump in the superfluid density occurs at the critical 
temperature Tkt [HE]- 
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Tkt for Superfluids 

By combining Eqns. ()B.lflj) and ()B.21|1 we are able to determine the critical 
temperature, T^t, for superfluids: 

1 2Kb '2Kb 



Tkt tt J tt hT 



(B.22) 



where the final form, in terms of the helicity modulus, is the most useful 
form. This can be written in lattice field theory (using lattice units Kb = 1, 
h = 1 and Lt = l/T) as: 



Appendix C 
Changes to QCD 

C.l Lattice QCD at ^ 

Lattice calculations cannot be extended to yU 7^ due to technical prob- 
lems. A quantum system, recast in Euclidean space-time, can be represented 
statistically in terms of the partition function: 

Z = j dUdi)di)e-^^^^'f'^^\ (C.l) 

where S = J2x ^QCD is the QCD action. Then, in order to calculate this, or 
related expectation values, 

< C >= ^ j dUdi)d^Oe-^^^'^'^\ (C.2) 

one can explicitly integrate out the bilinear fermionic dependence leaving 

Z = j dUdeiMe-^'^^^ (C.3) 

where Sg is the gauge contribution to the action and M = {Ip + mo) is the 
fermion kinetic matrix. Now, from the direct correspondence between Eu- 
clidean quantum field theory and statistical mechanics, it is possible to carry 
out the integral over gauge field configurations using Monte Carlo methods, 
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in which the highly peaked nature of e""^ is exploited by using it as a sam- 
pling weight. This is possible because Ip is an anti-hermitian operator which 
obeys chiral symmetry, so we have ,75] =0 and its eigenvalues come in 
complex conjugate pairs, ±iX. Therefore, the determinant of Ip can be writ- 
ten as a product of the eigenvalues, which in this case is real. Hence, we can 
use: det Jpe^^^ as the full functional weight. This is still a valid argument 
even with the inclusion of the chiral symmetry breaking mass term, as the 
eigenvalues of + mo are mo ± iA, which remain complex conjugate pairs 

However, when a chemical potential is introduced we have, M ^ Jp + 
f^o + f^lo- As the chemical potential term is hermitian, M then has com- 
plex eigenvalues^ making its determinant complex; and so the importance 
sampling weight becomes |detM|e*^e-^«. This can be shown as follows : 



i^l5{-ludu - /i7o + m)75?/^ 



(C.4) 



Therefore, 



75Mt(-;i)75 



det M{fi) 



det75Mt(-/i)75 



det M"f(-/i) 



= det M*(-/i) 



(C.5) 



^The determinant of Hermitian + AntiHermitian matrices is co: 
^ Using 7^ = 7^ and ^ -dl 



implex. 
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Then, for /i = 0, det M=det M* implies det M is real, but if /i 7^ 0, we have 
detM(yu)=det M*(— /i), implying det M is complex. 

Therefore, the case for /i 7^ leads to a complex phase, which in turn 
leads to configurations with large e~^^ cancelling and the effectiveness of the 
importance sampling weight becoming suppressed jTHI. 



In the NJL model, it can be shown that even with the introduction of a 
chemical potential term, the importance sampling weight remains real. Using 
the NJL Lagrangian, Eqn. ()8.7|) . we have: 



Now, we have^: 

^M(/i)^= (C75®r2)(C-S5®r2)(7^9^ + /i7o + m 



+ a + ij^f ■ 7i)ip 

= ip (C75(g)r2)C~S5(7i^5i/ + /^7o + "^ 

+ cr - ^757^ ■ 7f)r2^/' 

= ip (C75 O r2)C~^(-7^a^ - /i7o + m 

+ cr - i'y^T* ■ 7f)75 (g) T21/J 

= ^ (C75 ® T2){'^ldy + fi-fQ+m 

+ a - i'^lf* ■ Ti){C~^'^^® T2)ip 

= ^ (C75®r2)M*(/i)(C^S5®r2)^/' (C.7) 



Here we have det M = det M*, implying det M is real. 

^Using Cj^C^^ = —7* where the charge conjugation operator satisfies: C^^ ~ — C, 
and T2TT2 = — T* 



C.1.1 ;U 7^ in the NJL Model 



M ^ydy + yU7o + m + (a + i^^T ■ vf) 



(C.6) 
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C.2 Staggered Fermions 

The Euclidean action for a free fermion field ip, in 2 + Id, is given by 

S = d^xipip +7710)11;, 
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which, when discretised becomes 



2a 



Fourier transforming this to momentum space this becomes 



S = a-' 



(C. 



(C.9) 



(C.IO) 



(C.ll) 



(27r)3' 

a 

The inverse of M, where S = J ipMip, is found to be 

g^^^^ ^ --gY^ulv sin(ap^) +mo 

We can see that in the case of small momentum (long wavelength) limit, the 
small angle approximation implies sinapjy ~ apy and we recover the contin- 
uum Euclidean propagator. However, in general this represents 8 species of 
fermion. 

This "doubling problem" can be resolved by the staggered formulation. 
In this way we can reduce the number of species by a factor of 4, which 
can therefore by interpreted as the physical fiavours IHll- By making the 
transformation: 



^ iTiViTxx i^x ^ XxiTifiT 



(C.12) 



so that 7l)n,-iu7l)x±v vAx)xxXx±i>: where ?7^(x) 
action, Eqn. (jC.9|) . becomes: 



_l\xi+-+^^+\ then the 



El / \ iXx+u Xx-u) 
I ^''^^'^ 2a ^ ^oXxXx 



(C.13) 
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Now, M is in a reducible diagonal form, so disregarding all but one compo- 
nent of we are left with the reduced number of 2 quark species. 



Appendix D 

Hybrid Molecular Dynamics Code 



For the 2+ld Euclidean NJL model (in the Gor'kov basis), expectation values 
of an observable O as a function of the bosonic field configuration [$] are 
given by: 



where Z is the path integral (|8.24jl . As in statistical physics, such expec- 
tation values can be calculated using Monte Carlo or molecular dynamic 
techniques^ . In the most basic form we can consider that a statistical sample 
n = 1,2,..., A^} is used to approximate the full ensemble $; then the 
average of 0[$] is taken with respect to the randomly selected configura- 
tion of fields. However, as V det 2^[$]e~'^''°=[*] is highly peaked about certain 
configurations, it is then used as an importance sampling weight. 

In this work the expectation values were calculated using a hybrid molec- 
ular dynamics code 110], where the field configurations were evolved along 
classic trajectories with a time step of dt = 0.04. Then the expectation val- 
ues, (O), were conducted using complex Gaussian vectors (noisy estimators). 

^Where one evolves trajectories according to a stochastic process, the other evolves 
them along classic trajectories represented by Newton's laws 




(D.l) 
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This process was carried out for 300 trajectories. Then the final measurement 
was obtained by taking the arithmetic average of the various instantaneous 
values assumed by the HMD run. 

As an example we can look at the chiral condensate expectation value, 
Eqn. ()8.34|) . Once a trajectory is complete, and the end point is selected as 
the field configuration, the program calls the measure subroutine. In this 
part of the program a set of Gaussian vectors are defined as: 77^, where 
p = 1, ... ,4 are the isospin components and q = 1, . . . ,vol are the lattice 
sites; and if we average over the noise we obtain: 

Using this we can estimate a value for the chiral condensate. Starting with 
the matrix multiplication: 

V^A^'m(^__l^ '^^^rj = rj^Xrj (D.3) 

We can then average over the noise, as before, and we obtain: 

This will then be the instantaneous value assumed by the HMD run. Once 
this is carried out for the 300 trajectories, we then find their arithmetic av- 
erage and error. 

In this work we used an HMD code for the 3D four-fermi model with SU (2) x 
SU{2) symmetry, which was based on the algorithm in (41]. 



h6 



¥ 




(D.4) 
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